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^- _ Motivated by 3.9a evidence of a CP- violating phase beyond the standard model in the 

Qs ■ like-sign dimuon asymmetry reported by D0, we examine the potential for two Higgs dou 

in ' 

l/^ . blet models (2HDMs) to achieve successful electroweak baryogenesis (EWBG) while ex- 

m 



plaining the dimuon anomaly. Our emphasis is on the minimal flavour violating 2HDM, but 
our numerical scans of model parameter space include type I and type II models as special 
cases. We incorporate relevant particle physics constraints, including electroweak precision 



k> i data, b — )• sj, the neutron electric dipole moment, Rb, and perturbative coupling bounds 

C^ I to constrain the model. Surprisingly, we find that a large enough baryon asymmetry is only 

consistently achieved in a small subset of parameter space in 2HDMs, regardless of trying to 
simultaneously account for any B physics anomaly. There is some tension between simul- 
taneous explanation of the dimuon anomaly and baryogenesis, but using a Markov chain 
Monte Carlo we find several models within la of the central values. We point out short- 
comings with previous studies that reached different conclusions. The restricted parameter 
space that allows for EWBG makes this scenario highly predictive for collider searches. We 
discuss the most promising signatures to pursue at the LHC for EWBG-compatible models. 
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I. INTRODUCTION 

An experimental hint for beyond the Standard Model (SM) CP violation was recently reported 
by the D0 collaboration [lli|2|]. The like-sign dimuon asymmetry in semileptonic B decay has 
been observed to deviate from the SM prediction, initially with a statistical significance of 3.2 a. 
As this paper approached completion the statistical significance of this anomaly was reported to 
increase tli] to 3.9 a, while the central value shifted within the error band of the previous measure- 
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ment.^ This has motivated us to revisit electroweak baryogenesis (EWBG) in two Higgs doublet 
models.^ This measurement could indicate a new CP-violating phase contributing to the mixing 
of neutral B mesons or it could be a statistical fluctuation. Interpreting this observation as a sign 
of electroweak scale new physics is supported by a pattern of deviations measured from the SM 
in the B sector. "* This is also interesting as the measured CP violation of the SM is well known 
to be insufficient to generate the baryon asymmetry of the universe in electroweak baryogenesis 
(EWGB) scenarios. In this paper, we systematically reexamine the possibility of EWGB in two 
Higgs doublet models with new sources of CP violation in light of this experimental pattern. 

In any model with multiple scalar doublets, some symmetry (usually approximate) must be 
invoked to suppress flavour changing neutral currents (FCNC's). Frequently, treatments of multi- 
scalar doublet models impose ad-hoc discrete symmetries on the couplings of quarks and leptons to 



the two doublets, following Weinberg and Glashow lllOn . In this work, we focus on the framework 



of Minimal Flavour Violation (MFV) [1111 - 11311 . This approach allows a new physics sector with 
multiple scalar doublets that is naturally consistent with flavour constraints without such discrete 
symmetries. For the purposes of testing electroweak baryogenesis, the MFV 2HDM is sufficiently 
general that previously considered "type I" and "type 11" models can be thought of as special 
cases of the MFV Yukawa couplings. Therefore our analysis applies to a broad class of 2HDMs, 
even though our emphasis is on the MFV framework. Ref. lll4ll has recently considered EWBG in a 
2HDM which does not fit into any of these categories, but instead has couplings to the extra singlet 
that are engineered to reproduce the D^ anomaly while still respecting constraints on FCNCs. 



Electroweak baryogenesis in 2HDMs has been discussed in a number of previous papers [|l5- 
\2m and is the prototypical model for EWBG, where the dynamical source of CP violation is the 
simplest: it is the spatially varying phase 9(z) of the top quark mass."^ With the exception of ref. 



^ After this paper appeared, two new measurements were published. A new result from D^yft is still consistent 
with the dimuon anomaly, while a new result from LHCb is in tension with the dimuon anomaly. We discuss these 
results in more detail in Section IVTIl 

^ Only one linear combination of fields plays the role of the Higgs boson, so these should more properly be called 
"two scalar doublet models," but here we adhere to the customary convention. 

^ The measurements of Bg — > J/ip </> and B^ -^ tv also hint at the possibility of a new phase in i?s,d mixing (for 
a recent global fit and discussion see [4]). These deviations have the correct correlations to be part of a consistent 
pattern pointing to a new CP-violating phase. The deviation of B~ -^ tv has the most statistical significance — it 
is a 2.6 (7 deviation between the SM expectation |J4|] to the averaged measurements of B~ -^ tv performed at Belle 
and Babar 15|-|81]. See 191] for a recent discussion of B~ -^ t vin this model framework. 

■* Even though the MSSM has two Higgs doublets, it provides no such phase for rat at tree level because of the 
restricted form of the couplings in the Higgs potential. 



[12011 . these papers were written before there was a consensus on the proper way to compute the 
source term that appears in the Boltzmann equations needed to solve for the baryon asymmetry. 
Thus ref. [|20D gives the most complete treatment to date; yet it leaves much room for the im- 
provements that we undertake in the present work. We point out a heretofore unnoticed effect 
that systematically suppresses the magnitude of the phase 9{z), making it challenging to get a 
large enough baryon asymmetry. We note that actually solving for the bubble wall profile rather 



than merely parametrizing it, as was done in ll20ll . typically leads to smaller results for the baryon 
asymmetry. 

We also find that the particle physics constraints imposed here, but not in previous works, 
significantly reduce the parameter space that can succeed for baryogenesis. Using Monte Carlo 
methods we scan over the full allowed parameter space of the model consistent with these con- 
straints, rather than restricting it in an ad-hoc manner, as has been done in previous studies. 

The present work should be regarded first as a fairly general analysis of EWBG in 2HDMs, 
since as we will argue, the MFV framework is broad enough to also encompass some popular 
2HDMs (Type I and Type II) that have been considered before. Beyond this, we also consider 
whether it is possible to simultaneously account for EWBG and the level of new CP violation in 
the B sector suggested by the D0 anomaly. It will be shown that the two effects are largely un- 
correlated, due to the fact that we do not find significant EWBG associated with new CP violation 
coupling to b quarks; rather it must exsit in the top quark Yukawa coupling for this purpose. In 
fact, we will find some tension between simultaneously explaining the D^ dimuon observation 
and getting sufficient baryogenesis. 

The paper is organized as follows. In section HI] we describe the 2HDM model with MFV and 
fitting to the dimuon anomaly. Section Hill summarizes the various phenomenological and theoret- 
ical consistency constraints that we impose when scanning the parameter space of the model. In 
section|IV]we construct the one-loop finite-temperature effective potential that is used to determine 
properties of the electroweak phase transition. Section |V] describes the results of scanning over 
the model parameter space, using a Markov chain Monte Carlo (MCMC), to find models in which 
the electroweak phase transition (EWPT) is strong enough for EWBG. In section |VI] we study 
the baryon asymmetry generated by these models, finding that it tends to be somewhat below the 
observed value when the dimuon anomaly is reproduced, although exceptions can be found. Even 
when the dimuon anomaly is neglected, only relatively rare examples exist that can produce the 
observed baryon asymmetry. We then discuss in section IVIll the experimental prospects at LHC for 



confirming tlie dimuon anomaly and the MFV two Higgs doublet model, assuming that successful 
EWBG occurs. Conclusions are given in section IVIIIl 

II. MODEL 

We begin with some remarks motivating the choice of MFV as our principle for suppress- 
ing FCNCs. The mass scale suppressing the operators of interest enhancing Bs^d mixing (for 



perturbative couplings) is a few hundred GeV [12 ll . l22h . Models that seek to explain the D^ 



measurement with such a mass scale are strongly constrained by flavour changing measure- 
ments that agree well with the SM. New physics (NP) models with MFV are naturally consistent 



with such flavour constraints (in the Glashow- Weinberg sense [110(1) as the quark flavour group 
Gp = SU(3)uh X SU(3)dh X SU(3)ql is only broken by the SM quark Yukawa couplings gu/u 
defined as 

Cy = 9u ) u-M H^ e Qi - ^ ;. d.n H^ Q{ + h.c. (1) 

where 

1 

(2) 

-1 



2511 and can 



MFV models do allow new CP- violating phases that can lead to the D^ signal [|23 - 
possibly supply the required extra CP violation for EWGB. They provide an interesting framework 
for examining the relation between these experimental anomalies and the possibility of EWGB in 
SM extensions. It is also interesting to assess the feasibility of EWBG in multi-scalar doublet 
models with MFV, regardless of whether the dimuon anomaly is confirmed by future experiments. 
We denote by H the doublet that gets a vacuum expectation value (by definition the Higgs) and 
by S the doublet that does not. One can always rotate to a field basis where this is the case in the 
models we consider. Our discussion of the model will largely follow 125( 1 • The Lagrangian in the 
Yukawa sector is 

^Y = 9uiu,RH^ eQ{- gijid^nH^ Qi + Yu)u,nS'^ eQ{-Yu)d,RS^ Qi + h.c. (3) 

where flavour indices i,j are shown and summed over and color and SU(2)l indices have been 
suppressed. MFV asserts that any NP also has Gp only broken by insertions proportional to 
Yukawa matrices, so that Yjj-'^, Yj~,\ are proportional to g^\, g^^-. One can construct allowed NP 



terms by treating the Yukawa matrices as spurion fields that transform under flavour rotations as 



9u -^ Vu gu V^, go -^ Vd go V^^, 



(4) 



where Vu is an element of SU(3)ur, Vd is an element of SU(3)dr, and Vq is an element of 
SU(3)ql, i.e., the Yukawa matrices transform as gu ~ (3, 1, 3) and go ~ (1, 3, 3) under Gp. 
MFV can be formulated up to linear order in top Yukawa insertions, or extended to a nonlinear 



representation of the symmetry |26Ll27|] . For enhanced CP violation in Bq mixing we are interested 
in (at least) the second order terms in the expansion of the top Yukawa in MFV. It is sufficient in 
our initial discussion to only expand to next order in insertions of gu so that 

Vu gu\ + Vu dukiiaD^ i9uy^] + ■ ■ ■ , 



Y ^ 

^U i 



^ D i 



VDgDi + V'D9Dk[Wu)\ (9^ 



u) 



+ 



We decompose the second scalar doublet as S'^ = {S~^, 5*°), where 5*° 
scalar potential is 



(5) 

[SR + tsi)/V2. The 



V = - (h^^h,-- 

4 V 2 



m 



I {S^' S,) 



[m: 



H^'S, + h.c. 



+ Ai (H^'Hi) (S^'Sj), +A2 (H^' H,) {S^^ Si) + [XiH^' H^' Si Sj + h.c] , 

+ [X^H^' S^^ S, S, + X^S^' H^' Hi H, + h.c] + X(,{S^'Si)\ (6) 

where z, j are SU(2) indices. Here v ~ 246 GeV is the vacuum expectation value (VEV) of the 
Higgs field. Since we adopt the convention that the doublet S does not get a VEV the parameters 
m\ and A5 are related by, 



m; + Xt— = 0. 



(7) 



The spectrum of neutral real scalar fields consists of the Higgs scalar h = -\/25R(iJ°) another scalar 
field sji = y/2^{S^) and a pseudoscalar si = -\/2$j(S'°). However, these are not mass eigenstates; 
in the {h, sr, sj) basis the neutral mass squared matrix A^^ is 



^ ml X^v^ Xiv''^ 



M' 



X§v^ 



m 



H 



m. 



(8) 



where^ 



m\ = Xv'^/2 , mil = m| + X^v^ 



and 



^\ = ^'s- -^3^^ 



(9) 



We make A3 real by a phase rotation of S with respect H. We also define A4 = A^*^ + iX'^ and A5 = A|^ + lAg 
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with m| = ml + (Ai + X'^v^ jl. Note that rriu, mA is associated with s/j, sj. The mass eigenstate 
field basis is denoted as h', s'^, s'j and can be expanded in terms of the original field basis as 

h' = h- es^ sr - esj si, s'j^ = sr + e^^ h, and s'j = sj + esjh, (10) 

where we defined the expansion parameters 



^Sr = 2 2 ^^^ ^Sl 



rriH -mi m\- m^ 



NP 



The general Wilson coefficient for the relevant operator in the effective Hamiltonian H 

{y^^ Vthf C^^{mt) b'^ ql h^, q1 for perturbative Af '^ < 1 to leading nontrivial order in the MFV 



11 "*" r> 2 2 



c'^'im) = WoYy ' ' ^" ;"" " + ^%2 2 ' (11) 



where y^ = ^/2mt/v. In the above equation, the bottom quark mass m^ ~ 2.93 GeV is evaluated 
at the top quark mass scale. Here the Wilson coefficient was obtained under the assumptions 
rriH, rnA > m^, es^, £si ^ 1- In our convention the CKM factors are pulled out of the Wilson 
coefficient and appear explicitly in the effective Hamiltonian. Note that the Wilson coefficient is 
proportional to (?7^)^, not |?7^p. In general this parameter is complex, and the dependence on 
{tj'jjY introduces an extra phase into the mixing of Bg mesons through this Wilson coefficient. 

For EWBG, we will make use of the different but analogous phase that comes from the MFV 
expansion of the Yukawa couplings in the top sector. However to explain both of these anomalies 
simultaneously a number of phenomenological constraints common to the physics of both effects 
must be accommodated. For example, the mass scales of the new scalar states are relevant both 
for achieving successful baryogenesis and for explaining the dimuon anomaly. 

A. CP violation 

Since new CP violation is a central ingredient to this work, let us summarize the new CP- 
violating phases in our model. The new Yukawa couplings to the extra Higgs field S are beyond 
the standard model, and their phases cannot be removed by field redefinitions, which have already 
been used to push the phases in the SM-like Yukawa couplings (to the Higgs H) into the CKM 
matrix. For our purposes, the relevant new couplings are to the heaviest quarks, t and b, the first 
of which is important for baryogenesis and the second for the dimuon anomaly. 



Beyond these two new phases, the Higgs potential has four complex couplings, m^, A3, A4 and 
A5. Because of our convention that S does not get a VEV at zero temperature, eq. (|7]), m^ and 
A5 are linearly related, removing one of these phases. Moreover we can do a phase rotation on 
S (relative to H) to remove an additional phase; we choose to make A3 real. Thus there are two 
unremovable phases in the Higgs potential, in addition to the two new relevant phases in the 5* 
Yukawa couplings. We will show (see sect. IV Bl) that nonzero A5 suppresses the strength of the 
electroweak phase transition, so its phase will not play an important role for us. The effect of the 
phase of A4 is to allow for S to be complex inside the bubble walls of the EWPT, inducing a phase 
difference (psh between the neutral components of S and H. It is the sum of (psh and the phase 
of the new top quark Yukawa coupling, denoted by y?^, that will appear in the spatially-dependent 
CP- violating phase in the bubble wall, 9(z); see eq. (|6TI) . 

B. Relations to other Two Higgs Doublet models 

We will show that only relatively rare examples exist leading to a large enough baryon asymme- 
try in this model. We believe that paucity of models producing a large enough baryon asymmetry 
also holds for traditional type I or II two-Higgs doublet models as we now explain. Consider a 
type I or II two-Higgs doublet model where the matrix Yukawa couplings to Hi and H2 satisfy 

9h=9l = 0, type I, 

9u =9d = 0, type II, 
and gljjj denotes the diagonal Yukawa couplings of the scalar fields to up or down quarks respec- 
tively. Each of these Yukawa couplings is multiplied by 77^ and r]}j. In this notation, both fields 
generically get VEV's in the ratio {H2) /{Hi) = tan/3. The form of the couplings (fT2l) can be 
insured by Z2 symmetries, but in order to get a CP- violating phase in the Higgs sector for EWBG, 
these symmetries need to be broken; no unbroken symmetry distinguishes the scalar fields. 

If the symmetry were exact, then its action on the fields could be obscured by rotating to a basis 
in the Hi, H2 space. But since the symmetry is broken, nothing inhibits us from redefining fields 
to go to a new basis where Hi is along the symmetry breaking direction and H2 gets no VEV. By 
transforming the Yukawa couplings accordingly, we immediately find that the broken type I model 
corresponds to the general MFV model with the condition r]^/^ = —Vu/d ^'^^ 1^- Similarly, the 
broken type II model corresponds to the condition 77^ = —77^ cot (3 and rfj-) = —7]}j cot /3. We 
have set all other parameters v'u^Vu^v'd^Vd + ■ ■ ■ for both scalar fields to zero in this discussion 
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for simplicity. Our MCMC exploration of the parameter space would tend toward these relations 
between the couplings if they were favorable toward EWBG. The fact that they do not shows that 
the flavour structure of this model is not responsible for the difficulty of achieving EWBG in our 
analysis. 

This demonstrates that the MFV 2HDM is a general framework that includes type I, II models 
(in which the discrete symmetry is softly broken, as is always required to get EWBG) where both 
scalar fields couple to the quarks as special cases. For the purposes of only achieving EWBG, in 
fact there is no interesting distinction between the models, since as we will confirm the top quark 
Yukawa couplings by themselves give the dominant contribution to the baryon asymmetry. Thus, 
if setting aside the question of B physics, we can set 7]}^ to zero. Then the relations between tjIj to 
reproduce type I or II Yukawa structure can always be satisfied for some choice of tan /3. To the 
extent that the top quark source dominates EWBG, there is no significant difference between type 
I, II or MVF models. 

C. Fitting to the D0 and other B sector anomalies 

When attempting to link EWBG to the recent dimuon anomaly, we use the fit of i4\ to determine 
the new contribution to B^-Bq mixing (here q = s, d). The D0 result that deviates from the SM 
prediction at 3.2a {algi) and the SM prediction (A^^) are given by 

= -(9.57 ±2.51 ±1.46) X 10-^ (13) 

4l = (-3.10+°:^^) X 10-^ (14) 

The number of X bb — t- jX^ h^Y events is given by N^^ for example. The quoted a^^ is a 
combination of the asymmetry in each Bq, denoted a^j^. Each of these contributions to a^^ can be 
expressed in terms of the mass and width differences (M12, ri2) of the Bq meson eigenstates and 
the CP phase difference between these quantities (pq as 

IF'' I 

bq \^ 12\ ■ 1 /1 C\ 

«5L= mTT^ sm(/.g. (15) 

The model we discuss involves a NP contribution that includes a new CP- violating phase to M^2- 
It also allows other new CP- violating phases, all of which are interesting for our study, due to their 



11 



potential to drive EWBG. The effect of this model on Bg and B^ mass mixing is encoded in two 
real parameters, hg > and a^, by writing 

M!, = {Mf.f'' + (Mf^)^^ , {Mf.f = {Mf.f'' /..e^-^. (16) 

Frequently in NP models that obey MFV one has a flavour universal phase, so that hg = hd and 
o"s = cTd, and the best fit values are hq = 0.255 and 2aq = 180° + 63.4°. This scenario is argued to 
be a better fit to the current data then the SM in [|4|], which is disfavoured with a p-value of 3.1a. 
A recent update of this measurement by D^ with 9 fb^^ of data finds that the significance of the 
deviation has increased to 3.9a. We continue to utilize the global fit results when attempting to 
accommodate the dimuon anomaly and EWGB, as a new global best fit value is not available. The 
updated result is within the error band of the previous measurement that is incorporated into the 
global fit, so we expect our conclusions to hold for an update to the global fit incorporating the 
more statistically significant deviation. 

We treat perturbative QCD in the leading logarithmic approximation and evaluate the needed 
matrix elements of four quark operators using the vacuum insertion approximation at the b quark 
mass scale. The effective Hamiltonian is H^^ ~ (V^* VtbY C^^^rrit) b% q1 6^ gf (see Eqn. (fTT|) ) 
where a, /3 are colour indices; one finds 



{Ml,r ^ iy:^ \4,)^C^^(m,) ( -^ ) r^'flms,. (17) 



5 
24 

with the Wilson coefficient defined as above. Using the results of {2^ we have 

where r] ~ 0.84 is a QCD correction factor, and rj' ~ 1.45. The SM contribution is 

^C-K,^gM^s(|-). (19) 

where S{ml/M^r) ~ 2.35 [28]. Using these results we can scan over the allowed parameter 
space in the potential parameters Aj, couplings rji, and the masses of the new scalars, to search for 
examples consistent with successful EWBG while fitting the dimuon anomaly. 



D. Higher order terms in the MFV expansion 

Flavor breaking in MFV is based on an expansion in the insertions of the spurions Y^Yu 
and Y^Y^). These spurions can be inserted between any contraction of the Ql flavour indices. 
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Each term in the expansion where a spurion is inserted comes with an unknown parameter 
Vu^Vu^Vu ' ' ' ■ If these parameters are 0(1), it is not appropriate to expand to only leading or- 
der when considering the effects of the to p Y ukawa.^ Formally, one can treat these higher order 



terms in the elegant GMFV approach of ||27|]. In this section, we discus higher order terms in 



this expansion and their impact on the MFV 2HDM model. We adopt a less elegant, but simpler 



approach than [12711 for our constraints. 

We define new couplings, Cc/, C'u^ Co, C'd^ that resum parts of the series in the original expansion. 
The couplings C(7, Co are defined as the sum of all terms in the usual MFV expansion that do not 
involve flavour change, while the couplings ([j, ('^ are defined as the sum of all terms in the 
expansion in r]i (\^mt/vy that do lead to a particular flavour change. For couplings to neutral 
scalar fields this t ^ j flavour change is from one charge +2/3 quark species to another. For 
couplings to the charged scalar fields, the flavour change is from an up type quark to a down type 
quark. In each case, in MFV, the transition is accompanied by the appropriate CKM suppression 
of the standard model. Then for the neutral scalar fields one has for the top couplings 

Ct = yt{vu + vuyt+vuyt + ---), 

C = (n* H,) vt {v'u Vb+vuvt + ■■■), (20) 

where yt.fe = V^fntfi/v. For flavour change one must insert g\) go spurions as a correction to the 
gu coupling as flavour changing effects require the presence of both the up and down Yukawa 
matrices. These insertions, since they arise from a different spurion insertion, come with different 
parameters which we denote fj'^, fj'lj in the expansion of C,[^ . There are also g^^ g^ insertions in the 
expansion of (^^, but we neglect them in this expression. For the charged fields one has for the top 
couplings generating the t — )• 6 flavour change 

c;^ = iVu)yt{vu + Vuyt+Vuyf + ---)- (21) 

There are also glj go insertions in the expansion of Ci^ exactly as in the expansion of Cf- For the 
purposes of this discussion we neglect terms subleading in CKM or small Yukawa contributions 
to the effective couplings. The effective charged coupling multiplying the CKM insertion is the 



^ In the discussion above, we have expanded to the leading order required for the Bs^d mixing effect that we are 
interested in, avoiding these complications until this section in for the sake of clarity. 
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same as the effective coupling for the neutral field when no flavour change occurs. Furthermore 
Vtb ~ 1 to a good approximation. This identification allows one to directly relate the requirements 
of EWGB to the particle physics constraints that we impose. The procedure above can be repeated 
for the bottom quark couplings by defining (j^, C^^, Cb"^ Cb^- 

For the dimuon anomaly, we are interested in the effective flavour changing couplings involving 
the bottom quark for the new neutral scalar fields. We approximated this coupling with only the 
leading term in the previous section as {Vf^ Vts) rlj^ yb Vt- It is explicitly given by 

C = iV,l Vts)yb Wd yt+Vovt + ■■■)■ (22) 

so that the full result for the Wilson coefficient of interest for enhanced Bg mixing is given by 



c^^^K) 






2 



C. (23) 



iVtb^tsyb, 

where C is the quantity in curly brackets in Eqn. (fTT)) . The expression (^^^ appears in the precision 
measurement constraint formulae in the next section. When considering this scenario and fitting 
to the dimuon anomaly we will take 

C/iVbt) = fC/iV,tV,,) (24) 

with / = 0(1) as a generic relation between ^^°, which controls the strength of the dimuon 
anomaly, and (^j,^ which appears in the neutron EDM and 6 — )■ 57 constraints, as explained in 
sectionUnl The motivation for this choice is that these parameters get the same contributions from 
all the primed terms in the MFV expansion, while rjo only appears in (^^^. In the absence of fine- 
tuning (or some principle provided by the unknown UV completion of our theory), we expect this 
relation to hold with / ~ 1. However if rjo happens to partially cancel the other contributions 
to Cb^, or if the addition of all the spurion insertion terms sums to a suppressed value (compared 
to ?7£)), it is possible to get a relaxation of the loop-generated constraints relative to the _B-mixing 
contribution, hence smaller values of /. We will explore the dependence on / in the subsequent 
analysis. 

Aside from taking / ^ 1, another way of increasing the i?-mixing contribution without si- 
multaneously increasing the unwanted 6 — )• 57 and EDM contributions is to increase IA5I, which 
causes mixing between the scalar CP eigenstates. However this is not helpful in the current situa- 
tion. The problem is that the part of the function C that can be large enough to generate the dimuon 
anomaly (with smaller 77^) is proportional to (A^'^)^ and the first order phase transition required 
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,i?,/i 



for EWBG strongly prefers jAg ' | ^ 1. In physical terms, the mass mixing between H and 5* 
that is driven by A5 to allow the dimuon anomaly to be fit to (with smaller 77^) directly suppresses 
the first order phase transition required for EWBG. This is discussed in more detail in Section IVTl 



,RJ 



If IA5 ' I ^ 1, the mass and VEV eigenstates of the scalar fields approximately coincide, which 
makes discovery of the new scalar states challenging above the tt threshold as we will discuss in 
Section IVnl 



III. PHENOMENOLOGICAL CONSTRAINTS 

In this section we discuss the constraints on the 2HDM model that we impose when searching 
for parameter space with viable EWBG. The constraints are enumerated in Table I, which specifies 
the model parameters that are most directly affected by each one for the reader who wishes to skip 
the details. 



Constraint 


Constrained Parameters 


LEP/Tevatron Direct Search 


"ii,Ai,A2,Ct'^ 


RGE/Landau Pole/Unitarity 


A,mi,Ai,C° 


Neutron EDM 


Im[(C;^)*(C;^)^],m±(mi,Ai,2,3) 


6 — ;> S7 


Re[(c;^r(c;^)i,ic;^p,"^±("^i,Ai,2,3) 


Rb 


iCt'^nc;^p,m±(mi,Ai,2,3) 


EWPD 


"T-l, Ai^2,3 



TABLE I: Summary of constraints. 



A. Collider Mass bounds 



The kinematic direct production bound from LEP demands that the sum of new scalar and 
pseudoscalar Higgs masses obeys mn + rriA > 209 GeV. This is consistent with the assumption 
f^H, ITT' A > rrih made in deriving the Wilson coefficient in Eqn. ([TT]) . Moreover 



\{m]j 



m. 



m^ 



A1 + A2 



V > m 



h- 



(25) 
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The remaining direct production bound from LEP is on the charged scalar mass, m± = 
\/m{ + Ai/2 > 105 GeV. This bound is a purely kinematic constraint not dependent on tagging 
particular final states.^ 

For masses m± < rrit the constraints from the Tevatron can also be used. The CDF collab- 



oration 



m 



constrains Br (t — )■ S^ h) through subsequent decays of S^ — )■ cs, while the D^ 



collaboration 113 ih uses the subsequent decays S^ — > t^ u^.. The latter decay involves the cou- 
pling of the charged scalar to leptons, which is not directly related to the parameters of interest. 
Therefore we use the CDF result which for mf < 150 GeV gives the constraint 

Br {t -^ 5+ b) Br (5+ ^ cs) < 0.1. (26) 

For the 2HDM model with m± < 150 GeV we have Br {S^ — !■ c s) ~ 1 and this becomes 

8 TT miVf 



H '-t 



where r, ~ 1.3 GeV. 



B. Consistency Conditions on the Potential 

Although renormalizable, for our purposes this minimal model is best thought of as an effective 
low energy scalar sector with field content of a UV completion at a relatively low scale A ~ TeV. 
We will find that 0{\) parameters are required in the potential which can lead to a relatively low 
cutoff scale. For consistency we ensure that the couplings of the allowed parameter space do 
not approach a Landau pole or cause the potential to be unbounded from below up to A ~ 1- 
2 TeV ^ rrii, where rrii are the masses in the two scalar doublet model. We impose this constraint 
considering running and vacuum stability under the complete one-loop ROE flow of the couplings 
in the potential. 

The RGE running of the potential parameters can be derived directly by modifying the effective 



potential method of 113211 . For complex coefficients, within the convention that S is rotated so that 
A3 is real, the scaled beta functions are given in Appendix O Here the one-loop beta functions are 
defined in terms of the functions defined above as /3^. = (3^/16 tt^. In these expressions, g' is the 



^ The quoted bound of 78.6 GeV at 95%CL in 112911 is specific to the case Br {S+ ^ sc) + Br (S+ -^ r+ i^r) = 1- 



We impose the more conservative kinematic bound. 
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U(l) gauge coupling g is the SU(2) coupling of the standard model. The top Yukawa coupling of 
the Higgs field, and the complex Yukawa couplings of the new doublets are defined as 

l/* = V2^, Vt = C?, Vt = Cf (28) 

We also need the standard results of the one-loop running of the gauge couplings (above uit) and 
the Yukawa couplings; they are 

4, = 7g'\ 4 = -3g\ P, = -7gl ^ = y, ^- yl - ^^g'^ - ^-g' - 8gl 

rn n 17 n 1 

f^vt.t = Vt,b 



yt = yt 


"9 2 


- I9' - 8g! 



Q Q 17 

i^m ^i^\m -—9 --9 -8^3 . (29) 



Here we have neglected small mixing effects of the Yukawa coupling operators which we discuss 
below. For vacuum stability we assume that Ag > 0, which is sufficient for the field space direction 
along the S axis, considering the masses are constrained to be positive and large due to the bounds 
discussed above. For directions not along the h or (a particular) Si axis, it is difficult to analytically 
formulate the conditions for stability when A4 5 ^ 0, but it is easy to numerically check for runaway 
behavior, and this is the approach we take. 

From [133 we also have the unitary scattering constraint. As we scan parameter space and 
evolve the couplings of the model under the RGE we insist that the theory remain unitary as the 
couplings are evolved for mass scales up to the cutoff scale of the theory. However, we find the 
unitarity constraint has a negligible impact on the parameter space of interest; it is weaker than the 
demand to avoid Landau poles. 

C. Phenomenological constraints on the loop corrections due to charged scalars 

In this section we consider the constraints arising from virtual scalar exchange contributions to 
the neutron EDM, Z ^bb and b — )■ S7. In all of these processes, the dominant contribution comes 
from charged 5^ exchange because of the fact that the flavor change t ^ bin the loop is not CKM 
suppressed. 

1. Neutron EDM 

The effect of the exchange of virtual charged and neutral scalars on precision observables leads 
to significant constraints on models that would otherwise have given rise to successful EWBG. 
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FIG. 1 : One loop diagrams that induce an effective local operator of the chai^ged scalar field s^ to uji d^ 
from a coupling to ul d^. The largest contribution for the diagrams come from a top-bottom quark fermion 
loop in both cases. 



Some electric dipole moment (EDM) constraints on this scenario were discussed in [12411 . We 
focus on the EDM constraint that is not suppressed by small mixing angles or light quark masses 
that was discussed in BSll . For the neutron EDM we have 



Im 



2mtmi, 



h 



mr 



mr 



< 0.043, fg{x) 



xlog(x) x(x — 3) 



(a; -1)3 2 (a; -1)2 



(30) 



Naive Dimensional Analysis (NDA) is used to compute the neutron EDM matrix element in this 
estimate.^ To satisfy this constraint one generically has two choices: the parameter C^^ can be 
small, or the relative phases of (C*^)* and (C;,^)* can be tuned so that the constraint is satisfied.^ 
Either of these choices could conceivably be justified by further model building in the UV. 

Taking the parameter C,{^ small to satisfy the EDM bound generally requires |^^^| < 10"^ 
and the allowed C^^ decreases as C,^^ increases. One can also accomplish the suppression of C,^^ 
naturally by using MFV in a model where the new scalar doublet is not a flavour singlet; see 
[I25I1 . However, there is no exact symmetry in models of this form that distinguishes between the 



coupling of the charged scalars to ur di and ul dn; thus the effective parameters C^^ and C^^ are 
not independent and the corresponding operators mix. Therefore the radiative stability of such 
a choice is an interesting issue. An effective Ci^ is induced through the one-loop diagrams in 
Fig. [T] from the coupling of the charged scalars to ul dn proportional to (C;,^)*. Calculating these 



^ Although the uncertainty in the hadronic matrix element is significant and alternate estimates such as in II34I1 can 
weaken this constraint, we conservatively use the NDA estimate in restricting the parameter space. 

' In general (Ct )* and (C;, )* are the combination of many terms that can each individually have an independent 
phase in the MFV expansion. 
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diagrams one finds a contribution to the effective ur d^ Yukawa coupling such that 



5C^ ^ ^-^p^ f\x{l-x)\og 



iC^yy^ f\.n ^^,„J f^' 



\TX^ 







xvn^ + (1 — x)m^ 



+ ^fl/f I dx I \z(-\^{\^x^z)\oA^ 



1287r2 Jo Jo \ ' ^\A 



^(CfeJ^ /" rfx / \zil-x-z) 



\T\^ 



^0 



■'ASa^M,_U^ 



(31) 



Here A = — m^(x(l — x) — z^ — Axz — z) + ml{2xz + 2z'^ — z + x) — {1 — x~ z)My^r and 
we have retained the finite terms of these diagrams to illustrate the effect. For ^^^ > mt/nii, this 
contribution to the effective Yukawa coupling is large enough so that the induced effect on 6 — )■ s 7 
and the neutron EDM can be in conflict with experiment even if the coefficient of the operator 
Ufi di is chosen to be small at tree level. For this reason we restrict our choice of coupling to 
ICfe^l ~ 10 when scanning parameter space. 

For EWBG, it is advantageous to satisfy the EDM bound by tuning the relative phases. In this 
case Ct^ > C;,^ and one can have a strong coupling of the new Higgs S to the top quark. If on the 



other hand we were to suppress Cf" ~ 10 using the MFV model of |25fl or by choosing Q to be 



'± 



small, the small effective top quark coupling to the new field (and new CP violation) suppresses 
EWBG. 

2. Z — )■ 6 6 constraints 

Another important constraint comes from limits on the charged scalar masses and couplings due 
to modifications of the Zhh vertex. Z —^ bb provides direct bounds on the neutral and charged 
scalars.'*' The shift in Rb due to the virtual corrections of the new scalars comes about through the 
shift in gj^ = {gl^)sM + Sgj^, g^ = {gb)sM + ^Qh f^^ the b quark to the Z and we parameterize 
this effect on R^ by 

h )SM 

Jb )SM + \9b)sM 

We use the predicted SM values {J)^)sm = 0.0774, {gb)sM = -0.4208 which gives Rj, = 
0.21578 ± 0.00010. Using these results we have 

6Rb ^ -0.7S{6gt) + 0.U{6gl^). (33) 



XV r^oG {^ G\ f ^Sb i.9b )sM + ^Qb i.9b )am ^ 

dRb ~ 2 i?6 (1 - Rb) 1 —j—^ — . ,-B.2 I (32) 



'" See also II35L l36ll for recent discussions on Ri, constraints in models of this form focused on charged scalar ex- 
change. 
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Considering that the measured value given by the Particle Data Group is R^ = 0.21629 ± 0.00066 
B7ll we have the 1.3 a bounds -1.6 x 10"^ < 6 Ri, < 1.1 x 10^^ The shifts in the cou- 
plings, where we neglect suppressed corrections of the form rjiXj /IQn'^ due to mixing and 
ml/m\, Mz/m\ suppressed terms are given by 



Sg, 



-L 



IC 



'±12 



mf/m 



± 



Sg, 



T.R 



32 7r2 



'±12 
b I 

ic;^p 



m^/m\. 



rn^/m^ log (^mj/m\) 
{mf/m^. — 1)2 



5gt- 



(34) 



Here sw, cw are the sin and cosine of the weak mixing angle. We choose the renormalization 
scale yu = Mz- 



'± 



3. Constraints on Q through 6 — ;• s 7 



As explained above, we are interested in Q^ ~ 1 to allow for EWGB, and so we must in- 
clude terms in the h -^ sj constraint proportional to (Ct^)^- The constraint from BR(B — >■ 
s 7)i?^>i.6Gev at 95% CL ISSD is given by 



0.42 < -Re 



ich'iC^y I 



where from ll39h we have 



.(2) I ^ 



\\m?f^'' 



mr 



TTIj 



f\ i^) = ^7Z —^ + ^7Z TT^logW' 



x-1] 



4 x-1 



/f (^) = #^ + #^M^)- 



4 x-1 



2 x-1 



<0.12, 



(35) 



(36) 



D, Electroweak precision data 



As we scan parameter space we also restrict the mass spectrum of the scalars to be within the 
95% CL region of an electroweak precision data (EWPD) fit. The oblique parameters we use, due 
to the low mass scale involved, are the STUVWX parameters of ll40l l4in. The one loop corrections 



to the EWPD observables for the 2HDM are standard and are given in |25|] (with the appropriate 
rescaling). The fit we use is described in |42i]. We float the Higgs mass at one loop while testing 
the parameter set against the EWPD constraints and perform a joint fit to the Higgs mass and the 
new scalar spectrum. To do this we use the self-energy corrections due to the Higgs at one loop 
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given by 

,2N 91 



Ilzzip') 



IGtt 



^o^^/^) _ M^ i?o(/, m^ , M^) + i?22(/, mL M^) 



gf 



4 

Mml) 



IGtt^ cos^ 6*14/ 



Ml B,{p\ ml Ml) + 522 (/, m^, M|) 



(37) 



The one-loop functions used here are defined in ll42ll as are the STUVWX parameters 



This test is most sensitive to the mass splitting in the scalar spectrum. As the splitting violates 
custodial symmetry, it leads to important constraints on A3. For the MCMC scans we numerically 
interpolate from a grid scan of EWPD using the x^ measure defined through the cumulative dis- 
tribution function for a six-parameter fit. When we float the Higgs mass and scan the parameter 
space we interpolate between the discrete masses nih = {115, 130, 145, 160} GeV. For the new 
scalar masses we require niH, niA, rn^ < 700 GeV and vary the masses in this range over the 
allowed region considering all other consistency constraints on the scalar potential. The MCMC 
does not show a preference for masses exceeding this 700 GeV upper bound. 

Larger Higgs masses rrih are allowed in principle in these joint fits due to the positive AT 
contribution from the mass splitting in the scalar spectrum. We find that successful EWBG prefers 
a light Higgs. 

IV. FINITE TEMPERATURE EFFECTIVE POTENTIAL 

The properties of the electroweak phase transition and the bubble wall profiles needed for 
computing the baryon asymmetry require knowledge of the effective potential Vcs{H, S) at finite 
temperature. We compute it at one-loop order, along with the corresponding zero-temperature 
loop correction. The full potential can be expressed as 

Kff = ^Trce + V^W + V^CT + Vr (38) 

where Vrree is given by ^, Vr is the thermal contribution, Vcw is the one-loop Coleman- Weinberg 
(CW) potential, and Vct are counterterms, which for convenience can be chosen so as to maintain 
various tree-level relations for the scalar mass eigenvalues and VEVs. 
The CW potential is given by 

9i 4 f, ^1 3 



'^cw^Ei^'"? '"^-5 (39) 



647r2 ' \ Q^ 2 
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where gi is the muhiplicity of species i, rrii is its field dependent mass and + sign is for bosons and 
— sign for fermions. The field-dependent masses are given in Appendix|Bl Q is a renormalization 
scale, which we take to be v/^/2, where v = 246 GeV is the Higgs VEV. Some residual Q- 
dependence remains in our results since we do not perform a complete renormalization at 1-loop 
level including wave function renormalization. We include the Higgs bosons and electroweak 
gauge bosons in the sum, as well as the top quark. We work in Landau gauge where the Fadeev- 
Popov ghosts decouple at the one-loop level. 

The effective potential is a function of two complex fields, H^ and S^, but we have the freedom 
to remove the phase of H^ (the Goldstone boson mode) by making the appropriate SU(2) gauge 
transformation on both fields. Doing so simplifies our task by eliminating the unphysical degree 
of freedom. We thus regard Vcs as a function of the three fields h, s/?, s/ as defined below Eqn. (|7]): 

A. Counterterms and Goldstone bosons 

For convenience we introduce counterterms to preserve the tree-level relationships for masses 
and VEVs in the zero-temperature potential. In the {h, s^, sj} field basis we thus demand that 

dV/dh = dV/dsR = dV/dsj = at /i = f , sr = s/ = 0. Denoting the tree-level VEVs by {(pi), 
this requires 

(40) 



CT 



9Vcw_ 
d(t)i 



Similarly, to preserve the mass relations we require that 



d^Vc^ 



d(j)id(j)j 



9Vcw 



{<l>^) 



(41) 

{4>i) 



where (pi = h, sr, sj. There is a problem in principle with carrying out (HTj) . We work in the 
Landau gauge where ghosts decouple but Goldstone boson (GB) contributions are retained in the 
sum over species. In particular, intermediate GBs contribute to the mass of the Higgs, and formally 
their contribution to d'^Vcw / dh'^ is IR log divergent due to terms of the form {dml^/dhY In m^^,, 
where the prefactor does not vanish at the VEVs, but the argument of the log does. This shows that 
renormalizing the Higgs mass at zero external momentum with massless Goldstone modes is not a 
well defined procedure. Some authors choose to simply omit the GB contributions to the effective 
potential, but in the present study we find that this is not justified. Of course one can renormalize 
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the Higgs mass at any momentum scale p^ ^0 and a consistent procedure of renormalizing on- 
shell at p'^ = m\ was implemented in ref. ll43n . This cures the GB-problem because the external 
Higgs on-shell momentum flows through the loop of internal GBs contributing to the Higgs self- 
energy and cuts off the IR divergence that occurs when renormalizing at p^ = 0, Eqn. (|4T)) . The 



effect can be approximately captured by replacing m^^ — )■ mf^ in (141)) for the determination of 
the counterterms, where mf-^ is some IR cutoff. Since the IR divergence is only logarithmic, the 
effective potential is not greatly sensitive to the exact value. We find that mj^ = ml gives a good 
approximation to the more exact prescription of ref. ll43n . This is the procedure we adopt. 

In practice, not all possible counterterms are required. This is because we do not expect to be 
able to measure all of the couplings in the scalar potential in the near future, and therefore we 
do not need to relate all of them (such as \q) to observables. For our purposes, it is sufficient to 
include nine counterterms. 



VcT = Jq^^ h"^ ~ \5jj^ h^ + \5m\ (s^ + sf) + h {5m\j^ sr — 5m\j sj 
+ l6Xfh\sl - sj) - 5Xi h^ Sr Si + lh%5X^ sr + 6Xi sj) , 

which are fixed by the nine independent renormalization conditions set by eqs. 
formulas for the counterterm couplings defined in this way are given in Appendix lAl 



(42) 
Explicit 



B. Thermal corrections 



The unimproved pure thermal contribution to V^cfi is given by: 

dxx'^ln(lT exp(- ./5Tm2(^VS0)/T2 



Vt 






27r 



(43) 



where upper signs are for bosons and lower signs for fermions. Splitting Eqn. (l43l) into distinct 
sums over bosonic and fermionic species and expanding it to order 0{mf/T'^) one finds (see for 
example ref. MM)' 



mlT^ 
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771, 



B /j^'-B 



-Cb 



(44) 



24 127r 647r2 ^ T^ 

The lowest order potential (|43l) can be improved by resumming diagrams corresponding to inser 
tions of the thermal mass corrections. 



qO\ 



mi{H\ 5") ^ mi{H\ 5") + 5m. 



(45) 
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where Srn^ is of the form T^ times coupling constants. The Debye mass matrices can be obtained 
by computing the thermal self-energies in the high-temperature limit, but they can also be directly 
inferred from (|43l) . For a bosonic degree of freedom (pi the Debye mass matrix is 

^"^^3 2^ 24 90,90, ^ ^ 

k ■' 

where k runs over the relevant bosons and fermions in the theory, and m| are the field dependent 
masses given in Appendix |B}'. Fermions do not get a thermal correction. The thermally corrected 
mass mf + Smf in P31) denotes an eigenvalue of the full mass matrix m|^ + 5mfj. 

Strictly speaking the ring-improvement is a consistent correction only for the bosonic zero 



modes in the high temperature limit m/T <^ 1 ll44ll . correcting the masses only in the cubic terms 



in the bosonic expansion in (l44l) . However, we need the effective potential also for large m/T, and 
we need to be able to smoothly connect the two regimes. In the m/T >, 1 regime a cubic term with 
eigenvalues computed from high-T Debye mass corrections (|46|) becomes a poor approximation 
for the effective potential. The situation could be remedied by computing the Debye masses for 
arbitrary m./T, whence the Debye corretions would be exponentially suppressed at large m/T. 
This would be a lengthy computation but luckily a more convenient and mathematically equally 
consistent prescription is available: one implements (l45l) with the high-T expressions (|46|) directly 



in (|43] ) without making any high-temperature expansion ll45ll . Both approaches give essentially 
indistinguishable results in the high-T limit while the latter prescription also smoothly connects to 
the correct (vanishing) thermal correction at m/T ^ 1. 

To complete our approximation we also make the Debye correction within the 1-loop vacuum 
correction Vcw- Indeed, it is known that at 1-loop level the In(mf) terms cancel between the 
Vcw and the high-T expansion of Vt and we ought to preserve this property in our ring-improved 
potential. Log-terms corrected in this way do not induce spurious nonanalytic behaviour associated 
with the negative values of mf , which frequently occur for the Goldstone modes as well as for the 
physical Higgs in regions where the curvature of the potential is negative. Negative m^ induce an 
imaginary part of the effective potential due to the nonanalytic cubic T{m.'^y^'^ term, signalling 
an instability of homogenous zero-modes. Here we are only concerned with the real part of the 
potential obtained by replacing the cubic term with T|m^p/^. Numerically we use the high-T 
expansion where m./T is sufficiently small and smoothly match it onto a numerical fit to an exact 



" The sum over a given scalar of fermion representation can be written as J2k 9^'' 
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thermal integral for larger m/T. 

V. ELECTROWEAK PHASE TRANSITION 

We wish to identify models which have a first order electroweak phase transition, and which 
satisfy the constraint Vc/Tc > 1 for the Higgs VEV and critical temperature when bubbles of the 
broken phase nucleate. This constraint ensures that any baryons created during the phase transition 
are not significantly depleted by residual sphaleron interactions inside the bubbles. 

A. Finding Tc 

To determine the strength of the phase transition, we use the following algorithm. At suffi- 
ciently high temperatures one can usually find, through numerical minimization, a minimum of 
the potential that is in or close to the symmetric vacuum with vanishing Higgs VEVs. (We discuss 
exceptions to this statement below). At low temperatures one can of course find the broken mini- 
mum with (H^) = v/V2. The first step is to try to bracket the critical temperature T^ of the phase 
transition by subdividing an initial interval (Ti,T2) in the middle, Tm = |(Ti,T2). If only the 
symmetric minimum exists at Tm, then we replace Ti by Tm, whereas if only the broken minimum 
exists, we replace T2 by Tm- Continuing in this way, we will either reach an interval over which 
both local minima simultaneously exist, in which case we have bracketed Tc for a first order phase 
transition, or else the width of the interval becomes smaller than some cutoff (we take 0.1 GeV), 
and the transition is deemed to be second order. Once Tc is bracketed, one can accurately find Tc 
and Vc by minimizing the function /(T) = |Vsym — Vbr|, where V^sym and Vbr are respectively the 
values of the potential at the symmetric and broken minima. 

We find that for nonzero values of A4, the "symmetric" minimum in fact sometimes has a 
sizeable VEV Vs = A/2(|(if°)p + |(5'°)P). As long as Vg is small, sphalerons in the quasi- 
symmetric phase are still fast enough to induce baryogenesis. However if Vg is too large, the 
sphaleron interactions will be suppressed in front of the bubble wall and lead to a suppressed 
baryon asymmetry. We take the upper bound to be Vs/Tc < 0.4. This is the point at which the 
sphaleron rate in the "symmetric" phase starts to become significantly less than its value at Vg = 0, 
as can be seen from Eqn. (|67|) below. An example is provided by the parameter set nih = 122 
GeV, mi = 185.4 GeV, ma = 0, Ai = 1.335, A2 = 2.553, A3 = -1.291, A4 = -0.114 + 0.434z, 
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FIG. 2: Potential at T = Tc for a case that exhibits mild electroweak symmetry breaking in the quasi- 
symmetric minimum. The log of [V minus a constant near the minimum value] is plotted to exaggerate the 
barrier between the two minima. 

As = 0, Ae = 1, rjt = —0.096. The potential at the critical temperature is shown in Fig. [2l It has 
Vs/Tc = 0.34, and so still satisfies our criterion for sphaleron interactions being fast enough in the 
symmetric phase. 

B. Search strategy 



The scalar potential of the general two Higgs doublet model we investigate has many param- 
eters. While it is possible to scan them on a coarse grid to find many examples satisfying the 
particle physics contstraints discussed in section [nil very few of these models also have a strong 
enough phase transition, with Vc/Tc > 1. This is not surprising, since imposing constraints lowers 
the dimensionality of the allowed parameter space. For example, searching for points in an 12D 
space that happen to lie on a complicated 8D subsurface is bound to be inefficient if we restrict 
the search to points on a necessarily coarse grid in the full 12D space. Such a search could be 
made more efficient if we had an analytic formula for the constraint, but this is not possible for the 
EWPD. Table HI] shows the results of such grid searches for several values of nih. Later we must 
impose the additional requirement of getting a large enough baryon asymmetry, which is difficult 
to achieve. It is therefore desirable to maximize the number of models that satisfy the present 
constraints. 
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ruh (GeV) 


# models 


# strong EWPT 


115 


210,000 


92 


120 


195,000 


49 


130 


171,000 
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TABLE II: Results of grid searches of the 2HDM parameter space, showing the numbers of models that 
satisfy particle physics constraints, and of these the number having Vc/Tf. > 1. 

A more efficient search strategy is the Markov Chain Monte Carlo (MCMC) algorithm. Instead 
of blindly searching on a grid, this method favors points with a stronger phase transition, while 
still attempting to broadly explore the space. Starting at some point {Xj} in the parameter space, 
a chain of parameter sets is accumulated. A trial step is taken from the current point by varying 
the Aj randomly, with step sizes taken from a Gaussian distribution whose variance can be tuned 
for optimal performance. Let r be the ratio {vc/Tc)ncw / {vc/Tc)o\d, comparing the previous point 
to the new one. The new point is accepted into the chain if it has r > 1, and it is accepted with 
probability r if r < 1. In this way, a chain of models is obtained, whose probability distribution 
for Vc/Tc is proportional to Vc/Tc, thus favoring larger values of Vc/Tc. The occasional acceptance 
of points with smaller values of Vc/T^, is helpful to avoid getting stuck in a local extremum of the 
probability distribution, so that the parameter space is broadly explored. 

We modify the basic MCMC procedure by also rejecting the trial point if it has an unstable 
vacuum, or if it fails to satisfy any of the particle physics constraints described in section Hill We 
also adopt a trick to boost the acceptance rate: the trial point is alternately generated using the 
last point or the next-to-last point in the chain. Presumably the increase in acceptance results from 
the fact that the next-to-last point has already succeeded in leading to one new point in the chain. 
We set the Gaussian widths for the randomly generated steps to 10% of the current value for the 



parameters Ai, A2, A3 and ^J\m{\, and to 0.1 for the other dimensionless parameters Af , A4, Ae, 
C°, C° and the phase (p of C^- The width for steps in the mass of the Higgs is 5 GeV. 

The above description holds if one is only interested in generating models with large Vc/Tc. 
We would also like to favor models in which the baryon asymmetry can be large. Therefore we 
additionally compute the change of the CP-violating phase of the top quark mass across the bubble 
wall (777,^(2;) = \mt{z)\e^^^^'>), 

A^ = ^(^)r^ (47) 



27 

and test the quantity (A9){vc/Tc) instead of just Vc/T^. It then becomes consistent, within the 
philosophy of the MCMC approach, to impose Vc/T^ > 1 as another constraint to be satisfied by 
all models in the chain. Similarly, to search for models that lead to successful baryogenesis as well 
as the Dp dimuon anomaly, we take {A9)(vc/Tc)\Cb^ / (V^l Vts)\ as the quantity to be maximized 
by the MCMC. 

We have noticed that having nonzero IA5I tends to weaken the phase transition. (Recall that A5 
and m| are linked to each other by our choice of field basis in which (s) = at zero temperature, 
insured by Eqn. ([7]).) This can be understood analytically by considering a simplified version of 
the potential, 

V ~ Vh{h) + ^mls^ + \^sh{h^- v^) (48) 



Upon integrating out s, it becomes 



Vr^Vu{h)-^^h\h'-v'f (49) 



ni 



The second term always has positive curvature in the region h ~ v, tending to cancel out any 
barrier between the symmetric and broken phases that might arise due to finite temperature effects. 
We have compared MCMC runs in which A5 is fixed at zero to those where it is allowed to vary 
and find no significant difference with respect to the frequency of finding a given value of the 
baryon asymmetry. 

C. Results 

To illustrate some of the statistics of models with a strong phase transition, we consider a chain 
of 10^ models that pass all of the constraints. In this section and in section VI we use |?7^| as a more 
convenient parametrization for C,^ and similarly \riu\ for Q to avoid having to rescale by CKM 
factors and Yukawa couplings in the MCMC. (The relation between the rfq and C,q parameters was 
given in section HTPl ) We do not insist upon the dimuon anomaly here, but we allow |?7^| to vary 
so that we can check a posteriori what fraction of otherwise allowed models can be compatible 
with this potential constraint. We take the parameter / = 1 in eq. (I24l) so that the strength of the 
neutron EDM and 6 — )■ 57 are at a typical level given the assumed value of 77^. Fig. [3] (left panel) 
shows the distributions of ruh, the new scalar masses mn, ttia, rn± (where ttih/a were defined in 
Eqn. ^ and m^ are the charged Higgs masses), the dimensionless couplings IA4I, IA5I, Ag, |?7^| 
and \r]u\, and the quantities to be maximized (in addition to |?7^|): A6 (Eqn. (|47l) ) and Vc/Tc. One 
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FIG. 3: Left: distributions of parameters satisfying sphaleron and particle physics bounds, including ^^ 
and (^ , but not insisting on reproducing the magnitude of the observed D0 dimuon excess. mH,A denote 
masses of the new scalar and pseudoscalar Higgs bosons, respectively. Normalizations are arbitrary. Masses 
are in GeV. Right: distributions from MCMCs in which (^^^ = (^ ^ 0, and either omitting (heavy black 
bars) or applying (narrow red bars) the constraints from EWPD, b — ;• S7, neutron EDM, and perturbativity 
of couplings. Here Iwl) I^dI ^e proxies for IC?!; iCfe'^l- 

sees that relatively low Higgs masses nih < 140 GeV are necessary, the masses of the new scalars 
are less than 350—500 GeV for sr with mass Mh and the charged scalars (the pseudoscalar sj 
with mass Ma tends to be heavier, m^ < 500 GeV), and dimensionless couplings are typically 
less than 1 in magnitude with the exception of |?7^|, which we have deliberately pushed toward 
larger values in the MCMC. This parameter has no direct impact on the phase transition dynamics 
(since we do not include the b quark contributions to Vcs), but it does play a role in the Rb, b —^ 57, 
EDM and RGE constraints. 

To appreciate the impact of the new particle physics constraints we have imposed relative to 
previous studies, we also generated a chain of models that satisfy only the Vc/T^ > 1 requirement 
and the accelerator bounds on masses. The resulting parameter distributions are shown in the right 
panel of Fig. [31 They are much broader than those which incorporate the constraints. This explains 
in part why we find it more difficult to achieve enough baryogenesis in the 2HDM compared to 
earlier investigations. But there are other important reasons having to do with the production 
mechanism, as we describe in the next section. 
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VI. BARYOGENESIS 



To find the baryon asymmetry for a given model, once it has been established to give a first 
order phase transition and Tc has been determined, several steps must be taken. (1) The solution 
{h{z), sr{z), si{z)} for the fields in the bubble wall must be constructed. We will refer to it as 
the path taken by the wall in field space. (2) A set of coupled Boltzmann transport equations must 
be solved for the chemical potentials i-ii{z) and velocity perturbations Ui{z) of the particle species 
considered to be most relevant for creating a local CP asymmetry in the vicinity of the wall. (3) 
Some linear combination of the /ij gives the chemical potential for left-handed baryon number, 
/i^^, which biases sphaleron reactions in the symmetric phase in front of the wall. The integral of 
(xbl ill the symmetric phase region gives the baryon asymmetry, up to some proportionality factor. 

Our procedure is to carry out these steps for each model in the chains produced by the MCMC. 
We recall that these are models passing all the phenomenological constraints as well as the 
sphaleron bound Vc/Tc > 1, and also tending to have a large phase change A^, and optionally 
large values of ^j°. We weed out the models that necessarily give a small baryon asymmetry using 
a predictor. Statistically the dimensionless quantity 

g = — — -T (50) 

where L^, is the bubble wall thickness, tends to be proportional to the maximum value of the 
baryon asymmetry that can be produced. The correlation is shown in Fig. HI There are several 
reasons that the actual baryon asymmetry can fall below the maximum value predicted by q (as 
we will discuss below), but large values of ub almost always require that g ^ 0.15 (7^5/^5 obs)- 
Therefore if one wants to find the instances in a chain of 10,000 models that give the largest baryon 
asymmetry, it is not necessary to compute ub for all of them, but rather focus on those with the 
largest q values. Of course, q is much faster to compute than is hb- 

A. Bubble wall path and the CP-violating phase of the top quark 

In this section we discuss our procedure for finding the profiles of the fields h, sr, sj in the bub- 
ble walls, needed for computing baryon production. In principle, this task is complicated because 
the true bubble nucleation temperature Tnuc is somewhat lower than the critical temperature Tc. 
A complete account of the phase transition characteristics would involve finding Tnuc, computing 
the latent heat that drives the transition and evaluating the frictional forces exerted by the particles 
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FIG. 4: Scatter plot of baryon number in units of the observed density (n^/n^ obs) versus the predictor q, 
Eqn. dSO]). 

and fields on the expanding wall. Such a calculation would yield both the wall velocity and the 
Higgs field profiles over the transition region, which is what we need for a baryogenesis calcula- 
tion. This would be a difficult task however, and its results would still be uncertain because, e.g., 
reliable techniques to compute the interactions of the infrared gauge-field modes with the bubble 
wall do not exist. However, generically Tc — T^uc ^ Tc so that Vcs changes little between Tc and 
Tjiuc and one might expect that the true profiles are reasonably well approximated by the solution 
interpolating between the degenerate minima at T = T^. This is the approximation we shall adopt 
here, and the wall veolocity v^ is left as an external free parameter. The wall profiles are then 
found by numerically solving the equations of motion 



■'cff 



dz"^ 



(51) 



at the critical temperature, where z is the direction transverse to the wall, subject to the boundary 
conditions that the fields approach the broken or symmetric minima as z — )■ ±oo. 

1. Gauge dependence of bubble wall profiles 

The coupling between the neutral scalars and the Z-field poses a complication in solving for 
the bubble wall profiles. This is because the co variant derivative terms \Dz(t)i\^ provide a source 
for the classical Z-field, and so one is faced by the problem of minimizing the Hamiltonian: 



dz \D,Hy + \D,Sy + V,siH", 5^Te) + -M'^Z^Z'^ 



(52) 
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where D^ = dz — i{g/2cos9w)Zz and dots indicate other Z-dependent terms. Writing -\/2iJ° = 
fie^fh and v^S*" = se*"^^ and using the fact that the effective potential depends only on the phase 
difference (fh — ^s one can show that 

d, {f, - MlZ^) = , (53) 

where M| = g'^{h'^ + s^)/(4 cos 6'^) and j^ is the current sourcing the Z-field: 

f, ^ —^ {h'd^^n + s^d>^^s) . (54) 

2 con Oy^r ^ ' 

Since all classical fields should vanish far away from the wall, we find that 

while all other components vanish. Recall that in our formulation of the effective potential, we 
have chosen a gauge where the Goldstone boson vanishes, hence ^ph = 0. This leads to a non- 
vanishing source (|54|) for the Z-field since ^ps is generally nonzero. Of course Z-field is gauge- 



dependent, and by an alternative choice of the basis for phases [|46ll we could make j^ = 0, and 
in this gauge Z^ = 0, by Eqn. (l55l) . The price to pay for this simplicity is that one would have to 
keep track of four fields in the effective potential rather than three. 

2. Spatially varying phase of mt 

A key quantity for computing the baryon asymmetry in the current scenario is the spatially 
varying CP- violating phase of the top quark mass, given a priori by 

mAz) = ^e'^^ (h + ^ se^^^A , (56) 

V2 \ yt J 

where (psh = Ps — ^h- As we have observed, because of the U(l)y gauge symmetry, p^ and ps 
are not both physically meaningful quantities, so it is interesting to see how both gauges discussed 
above lead to the same phase for the top quark mass (especially since different authors have chosen 
different gauges in previous treatments). First, in the Z = gauge we can solve for the phase ph 
using the constraint j v = 0: 

dzPh = —nr, — ^dz'^sh ■ (57) 

Ph can be reconstructed after the fields /i, s and ^psh are determined using equations following from 
the reduced Hamiltonian: 

j dz {^-{dzhf + ]^{dzsf + \-^^{dzP^sh? + V,^{h,s,psh.TS\ , (58) 
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where V{h, s, cpsh) = Vcs{h, sji = s cosOsh, sj = s sinOsh, T^). (We can always shift cp^ into the 
phase of S in Vcs because it is invariant under global SU(2)xU(l)j^ transformations, and contains 
no derivatives of the fields.) 

Alternatively, if one sets (ph = 0, it is necessary to deal with the nonvanishing Z^-field: 

Z4^) = - ^^j^s'd^ips. (59) 

This field induces an additonal CP- violating force acting on the top quark due to the gauge interac- 
tion (g/2 cos 9w) i^lbt. This term can be removed by a local axial transformation t — )■ e*"^ "^^/^t, 
which however reintroduces an additional overall phase cpz into mj: 

mt(z) -^ ^ e'^^ (h + ^se*^M . (60) 

V2 \ yt J 



It is easy to see that this transformation removes the interaction term with Z^ given by Eqn. (1591) 
only if dz^fz = —jir^dz^s which is identical to Eqn. (l57l) in the gauge (ph = 0. Denoting 
Ct Ivt = ICt'/l/tk^'^'S we find that the total phase of nit is given by 

(I — |ssin(o9n + o9s/i) \ 
^^0 ' . (61) 

h+ |^|SC0S(V9^ + fsh) J 

It is noteworthy that both contributions to the overall phase are suppressed by powers of s/h: there 
is no phase if the bubble wall profile remains strictly along the h direction. Bending of the path in 
the s directions is necessary for generating any CP violation in the wall. 

Although it might appear from (l57l) that one more integration is required to solve for 9{z), 
in fact we never need 9{z) itself, but only 8^9, to evaluate the CP-violating source term of the 
Boltzmann equations, to be discussed below. Eqn. (|6TI) thus constitutes a sufficient solution for 
the phase of nit in the bubble wall. 

3. Numerical solution for bubble wall profiles 

The equations of motion can be solved numerically by shooting or by relaxation. However, for 
the purposes of doing an MCMC search, a simpler and faster approximate method is desirable. We 
first find the straight line in field space that connects the two VEVs and then locate, for each point 
on this line, the position along the plane perpendicular to the line where the potential is minimized. 
This deforms the line into a curve, which is the path of steepest descent through the saddle point. 
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FIG. 5: Left: the phase of S''^ -field computed from the full EOM (solid black) and from the minimization of 
V (red dashed) as a function of zT corresponding to model 1 in table |IVl Right: same for the zT derivative 
of the phase. 

One then solves the constrained one-dimensional problem along this path in the field space to 
determine the profile. This procedure is valid if the path does not curve very much (which we 
typically find to be true), but it tends to allow the path to curve more than the actual solution, since 
it ignores gradients of the fields when it determines the path. The real solution is systematically 
smoother than than the approximate one, which tends to overestimate the baryon asymmetry by 
a factor of <, 2. The approximate method is thus sufficient for searching good candidates in our 
MCMC-runs, but for accurate results we recompute profiles for all models that pass the constraints 
by solving the equations of motion by use of accurate relaxation methods. 

In Fig. [5] we show the solutions for the phase i^^^ (left panel) and it derivative (right panel) as 
a function of zTc using the minimization of the potential (dashed curves) and the full equations 
of motion (solid curves) for the second model listed in table |lVl As claimed, the true solution 
is smoother than the one determined by constrained minimization of the potential. In Fig. [6] we 
display the profiles \h, sr, sj} as a function of zTc for some other typical models. Based on 
these figures the wall thickness appears to be of order 10/T. This is an important parameter for 
electroweak baryogenesis and although we do not approximate the wall profiles by any analytic 
ansatz, it is convenient to characterize the wall thickness L^ by fitting the total VEV v{z) to the 
form if c tanh(z/L^). In this way we find that L^Tc ~ 6 — 12 for most models in our chains. The 
distribution of L^ values is shown in Fig. U\ 



34 




0.8- 
0.6- 
0.4- 
0.2- 



-0.2 



^0 



-20 




-10 10 

zx T 



20 



30 



FIG. 6: Examples of bubble wall profiles for the fields h, sr, sj. Left: a typical profile; right: a profile 
leading to an atypically large CP- violating phase 9{z), and using model parameters that are ruled out by Rb 
and 6 — )• S7 constraints. 

4. Comparison with previous results 

We have noted that the CP- violating phase 6{z) of the top quark mass is suppressed by the ratio 
s/h in the bubble wall, and that s/h is typically small for models that are not ruled out by particle 
physics constraints (see Fig. (6]). This is in contrast to the earlier analysis of ref. ll20ll . Contrary 
to this one and earlier papers, we are finding a systematic suppression in the magnitude of 6{z) 
that goes beyond the smallness of CP-violating phases in the potential. Namely, since {S) = 
at zero temperature, but 9(z) vanishes if S were to remain zero throughout the bubble wall, any 
nonvanishing 9{z) comes about as a finite-temperature effect, i.e., (S) being no longer zero in the 
broken minimum at finite T. On the other hand, ref. ll20n did not observe any such suppression. The 
reason is that ref. ll20ll did not actually solve for the bubble wall profiles, but made the assumption 
that the two Higgs fields remained equal to each other in the wall, \Hi\ = \H2\ (which would 
translate into l^l = in our field basis). This assumption is only true if Hi -H- H2is a symmetry 
of the Lagrangian, which ref. [|20|l claimed to be the case due to a choice of relations between the 
parameters of the Higgs potential. It would require A4 = A5 = in our language. But it is clear in 
that case that \S\ remains strictly zero in the bubble wall, so that 9(z) = 0. In the language of ref. 
[I20I1 . the CP-violating phases in the Higgs potential should vanish in order to have the symmetry 
ij"i -H- H2 leading to li^il = \H2\ in the wall, but this would also have made 9{z) = 0. In short, ref. 
[|20|] made inconsistent assumptions which obscured the suppression of 9{z) that we have made 
explicit. 
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FIG. 7: Distribution of bubble wall thicknesses L^, in units of T^ ^ 



B. Transport equations and source term 
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The next step in determining the baryon asymmetry is to solve for the chemical potentials of 
relevant particle species in the vicinity of the bubble wall, that are induced by the top quark mass 
phase 0(z). One approach to deriving these equations is to start with the Boltzmann equations 
and to perturb the particle distribution for species i away from its equilibrium form by allowing 
for a small chemical potential fii and a perturbation 5fi that describes the departure from kinetic 
equilibrium ll47n : 



h 



exp{-(3'y^{E + Vy,p^) - /ij) ± 1 



5f^ 



(62) 



where v^ is the bubble wall velocity and 7^ 



"^/^. Linearizing in /ij andwj = (5/jW^),the 



perturbation to the particle's fluid velocity, leads to a set of coupled first order transport equations 
for fii and Ui 147}]. This procedure has been carried out for the 2HDM in refs. I12d . 14811 and we adopt 



the same transport equations as in [|20ll which were tailored to two Higgs doublet models. Here the 
CP- violating phase of m^ provides the source term and the species that are followed are the quarks 
^L, tR., bL and the Higgs bosons, which are assumed to have a common chemical potential fXh and 
velocity perturbation u^. 

The transport equations depend upon the bubble wall velocity v^, which we take to be an un- 
determined parameter. In the computations that follow, we have chosen v^ = 0.1 for definiteness, 
but we have checked that the final results for the baryon asymmetry do not depend strongly on 
this choice over a range of reasonable values. In fact the most important t>^„ -dependence is that 
which multiplies the source term, which we discuss next. However, in the end this dependence 
gets canceled by an explicit factor of l/v^ in the relation between the chemical potentials and 
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the baryon asymmetry, eq. (1661) below. The residual dependence of the result on v^ comes only 
through its appearance in subleading terms in the transport equations, which explains the relatively 
weak effect. For example we find that the baryon asymmetry increases by only 20% when taking 
Vw —7- 0.05. (Ref. uM which uses a different form for the diffusion equations finds a stronger 
dependence.) 

1. The source term 

A crucial ingredient in the transport equations is the inhomogeneous source term, which after 
many years is still a matter of controversy in the EWBG community. The source term adopted 



by boh and by us encodes the fact that particles, and here the top quarks in particular, experience 
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M- The 



a CP-violating force while traversing the bubble wall, where 6{z) is varying [|47l . 
corresponding force term in the Boltzmann equations thus gives rise to the (scaled dimensionless) 
source 

St = v^ {K^{xt){xle')' - K^{xt)xlxl'e') , (63) 

where primes denote d^x and Jig g are dimensionless functions of a;^ = \'mt\/T arising from phase- 
space averaging of certain kinematic variables, ^^ plotted in Fig.[8l 

Since the typical wall thickness is significantly larger than T~^, the derivatives in St are a source 
of suppression of the baryon asymmetry. A source term with fewer derivatives could possibly give 
larger results. Such contributions arise from corrections due to CP-violating dispersion relations to 
collision integrals in the Boltzmann equations. For example the Higgs decay to tt pairs contributes 



a source 15311 : 



S^ = -v^^ U QiMH - 2\mt\)xle' , (64) 

where Mh is the mass of a Higgs boson, is the Heaviside function and /^ is a thermal kinematic 
function of {Mh/T, \mt\/T) typically of order ~ 0.1. However because of its small coefficient 
we find that the contribution of S^ to the baryon asymmetry is always at least 10 times smaller 
than that of 5*4. 

The controversy alluded to above stems from the fact that there are co mpeting formalisms for 



deriving source terms for the transport equations, that originate from refs. [1541 . l55h . These attempt 



to capture the quantum mechanical aspects of the particle interactions with the wall, rather than 



'^ Our Ks,9 ^e related to those of p20|] by simple T scalings that make them dimensionless. 
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FIG. 8: Left: thermal functions Kg g appearing in the top quark source term, Eqn. ( [63l ). and A'l which 
appears in the Bl chemical potential, Eqn. (1651 ). Right: an example (model 1 in table ITVl) of the source term 
as a function of zT. In this case the baryon asymmetry computed from full and ID approximation to the 
field equations agree to 1%. 



the semiclassical behavior described by the sources (I63H64I) . A difficulty with this approach is 
that rather drastic approximations need to be made in order obtain concrete, tractable expressions, 
notably expanding propagators to leading order in VEVs of the spatially varying Higgs fields (the 
"VEV insertion approximation"). The resulting source term has only a single derivative, unlike 
the semiclassical force source term we have adopted. It thus gives a larger estimate for the baryon 
asymmetry. However we do not believe that it represents a controlled approximation, so we adhere 
to the semiclassical approach. A more complete transport formalism, capable of accounting for 
the nonlocal coherence effects associated with the quantum interactions, is being developed in 
refs. pal . There has been recent progress in improving upon the VEV insertion approach in 



I157L 



580, although so far only with source involving scalar fields. 



2. Transport equations 



To solve the transport equations of ref. [|20h accurately, it is necessary to use the relaxation 
method, since the asymptotic behavior of the solutions in the large \z\ regions cannot be deter- 
mined analytically. On the other hand, shooting from large values of ±|2;| toward ^ = is compu- 
tationally much more efficient. The problem with shooting is that in principle one needs to know 
the analytic behavior of the asymptotic solutions in order to correctly specify the ratios /ij/wj. 
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FIG. 9: Typical example of solution to transport equations for the chemical potentials of ti, t\, bi and 
Higgs bosons. This example (model 1 from table Hvl) gives a baryon asymmetry which is 1.1 times that of 
the observed value. 

since only half of the independent variables can be input as free parameters (to be determined by 
smoothly matching the solutions at z = 0) at either boundary. For example, one could freely vary 
Ui at each boundary while holding /i^ fixed to find a solution, but one does not know beforehand 
what values of /ij to choose. If the starting points are sufficiently far from 2; = 0, then it might be 
a good approximation to take /ij = at the boundaries. But since the solutions are exponentially 
decaying functions of \z\, shooting is not numerically stable if the boundaries are taken to be so 
distant that setting /ij = would be a good approximation. 

Of course the extra computational burden of relaxation is not an issue if one only needs to solve 
the transport equations for a small number of models, but to check thousands of models generated 
by the MCMC, a compromise is needed. We find that shooting typically gives an estimate for the 
baryon asymmetry that is accurate to within a factor of a few. We reserve relaxation for refining 
such estimates for potentially interesting examples. Typical examples of solutions for iii{z)/Tc 
and Ui{z) from the relaxation method calculation are shown in Fig. |9l 

C. Baryon asymmetry results 



Once the chemical potentials around the bubble wall are known, they determine the overall 
asymmetry in left-handed baryon number, jib^ > which biases the sphaleron interactions. Following 
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[|20l1 . we take this to be 

fiB, = 1(1 + ^Ki^t)f^t + Kl + 4i^i,6)/"b - 2ifi,t/Utc (65) 

where Kn is a dimensionless function of mj/T shown in Fig.[8l (Notice that A' i f, = 1 since b is 
nearly massless.) The baryon-to-entropy ratio can then be calculated as 

405r,ph 



Vb 



y"rf^/iB,/sphe-^^^-'^l^l/(^''») (66) 



where Fsph ~ 10~^Tc is the rate of sphaleron interactions in the symmetric phase, the exponential 
factor takes care of washout in the symmetric phase in case v^ <^ 1, and f^ph takes into account 
the z-dependence of the actual sphaleron rate as it turns off in the broken phase. We estimate /sph 



as follows: in the broken phase where the Higgs VEV v = y 2(|iJ°p + \S^\'^) is nonzero, the 
sphaleron rate is given by'^ Fbp ~ A^e~^°''/^. We fix N by equating Fbp = H^^ at v/T^ = 1 
(this encodes the usual sphaleron wash-out constraint v/T^ > 1), giving Fbp = iJc~^e^°*^^"*'/^'=^ ~ 
2.4 Tc e''^'^'"^^", where we used H^ ~ lO^^^T^. We now assume that this equation holds also 
for V = v(z) within the wall, and define the overall sphaleron rate to be min(Fsph, Fbp(2)) = 
Fsph/sph(-2) which gives 

/spb(.)=minfl, ^e-^°^(^)/^') . (67) 

V J- sph / 

It is important to use a realistic estimate for fsph rather than a simple step function, because we 
find that the integral of hb^ over the symmetric phase can be much smaller than that of /i^^/sph- 
This method also treats correctly those cases (and automatically cuts out the most extreme) where 
the VEV is nonzero even in the symmetric phase. The effect of the fsph modulation is illustrated 
inFig.[IOl 

We have searched the parameter space for models likely to provide a large enough baryon 
asymmetry and, optionally, a large enough dimuon excess in addition. We consider two cases. In 
the first, we set (l^ = Q"^ = so that the constraints from Rh, b —^ S'j and the neutron EDM 
are relaxed as much as possible without having to fine-tune the parameter / in eq. (l24l) . In this 
case we have or course given up on trying to explain the dimuon excess. We find that it is difficult 
to achieve a large enough baryon asymmetry even without this extra demand. A chain of 10,000 
models designed to achieve large Vc/Tc and 9{z) yields no models that have r]B significantly greater 



'^ To be more precise the factor 40 in the exponent should be replaced by a function depending on the couplings in 
the scalar sector, but that would not change our results essentially. 
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FIG. 10: Example of z-dependence of fiB^, /^Bi,/sph and their integrals, showing the importance of the 
sphaleron rate modulation factor /gph- Symmetric phase is on the right. 

than the observed value, and only a handfuU which come close to the observed value. Rarely we 
find examples where the Landau pole occurs above the scale A = 3 TeV; working models with 
A = 1 TeV are much more frequent. The distribution of values of 77^ for these 10,000 models 
is shown in Fig. [HJ It is remarkable that the upper limit on the baryon asymmetry happens to 
be so close to the observed value. In the figure we also show the distribution of tjb in a chain 
where the constraints from EWPD, R^, b ^ s^ and the neutron EDM are omitted; in that case it 
is easier to reproduce the observed baryon asymmetry, although there is still relatively little room 
for producing too great an asymmetry. Parameters for a few working models are given in table 
Uni They share several properties: light SM-like Higgs masses ~ 117 — 125 GeV, new neutral and 
charged Higgs masses near 300 GeV, and values of \Cf\ ~ 0.5 that saturate the constraint coming 
from Rf), with C° having a large phase 9. 

In the second case, we bias the Monte Carlo to produce models with large values of |Ch°| to try 
to simultaneously explain the dimuon excess. In this scenario, while it is more difficult to achieve 
a large enough baryon asymmetry than in the previous one, the results are rather sensitive to the 
choice of / in (l24l) . with / — ;■ giving the most relaxed constraints. Similarly to the previous 
case, with a Landau pole cutoff of 2 TeV, we find no viable examples but at 1 TeV and we find 
several (out of iC models in the chain) that have a large enough baryon asymmetry. If / = 1, 
the actual values of 77^ ~ 1-2.5 fall below the central values 3.6-4.6 needed to agree with the 
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FIG. 11: Solid bars: frequency of tjb values from chain of 10^ models that satisfy all constraints. Open 
(red) bars: frequency of tjb when the Rb, 6 — )• S7 neutron EDM and perturbativity constraints ai^e omitted. 
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0.10 
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1.96 


0.52 
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0.52 


3.41 


1.07 
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1.20 
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2.26 


0.97 


-0.74 


0.12 


0.22 


2.14 


0.43 


3.68 


1.04 


1007 


1.18 



TABLE III: Parameters of sample models with Q^ = Cb — ^ (hence r]'^ = 0, and no dimuon anomaly) 
and sufficient baryogenesis, and a Landau pole at scale A. Masses/energies are in GeV. 

D^ dimuon observation. However a modest relaxation of the constraints by using / = 0.5 leads 
to more favorable results, where we find ~ 10 models with the desired baryon asymmetry and 
dimuon signal. (The distribution of models satisfying both constraints is displayed graphically in 
the final figure in section IVIIIl ) These examples are similar to those with (l^ = ^^^ ^ in that the 
SM-like Higgs must be light, while the new scalars exceed 300 GeV in mass. 



Vn. EXPERIMENTAL TESTABILITY 



We have shown that a direct link between the new phase hinted at in B physics data and EWBG 
cannot be made in 2HDM models in general. EWBG depends upon CP violation in the scalar 
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0.08 
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1.14 
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316 


2.36 


1.35 


-1.06 


0.052 


-0.215 


1.40 


0.381 
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1.24 
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1.04 
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0.046 


-0.168 


1.73 
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4.41 


3.7 
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TABLE IV: Sample of model parameters giving sufficiently lai^ge \r]'j^\ (hence ^['^) and sufficient baryogen- 
esis. ri'^0 is the value needed to match the central value of the dimuon anomaly. 

potential and top quark Yukawa couplings, while enhanced Bg mixing involves new couplings to 
the b quark and phenomenological constraints severely limit the latter couplings. Nevertheless, the 
prospects for direct collider confirmation (or refutation) of the new phase hinted at in the dimuon 
data and the new scalar degrees of freedom we have discussed are very strong. In this section we 
briefly discuss some of these prospects. 



A. Confirmation/refutation of the dimuon anomaly 

n 

The hint for new physics near the electroweak scale indicated by the D0 collaboration [20 that 
has at least one new source of CP violation, is consistent with other direct measurements that 
are sensitive to the same phase in Bg mixing and use flavour tagged B^ -^ J/\E' decays. With 
6.1 fb~^ of data, D^ [|59|] finds a hint of a new phase in decays of this form; their results show 
a larger deviation from the SM than the equivalent measurements by CDF [|60|]. Both results are 
statistically consistent. A more recent measurement that^has appeared since the original version of 
this work includes 8 fb"^ of data from D^ reported in [3]. A method for checking the consistency 



of these results can be found in [|21 



6III : one can compare the fit values of AF^ and Sf^^ from the 



two measurements. The relevant relation is 



|AF,| ^ - [0.20± O.llps-^] yjl - Sl^/S^^, 



using 5ms = 17.78 ± 0.12 ps"\ a|^ = -(4.7 ± 4.6) x 10"^ and a^^ 
0.093) X 10^2. Using S^^ = -0.55 ± 0.35 one finds that [3] implies 



C2 



-O.lliO.lOps^^ 



(68) 

-(0.787 ± 0.172 ± 

(69) 
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which is consistent with the updated value derived from the dimuon anomaly. It is encouraging that 
D^ has just announced that with 9 fb^^ of data the significance of this deviation in the dimuon 
measurement is now 3.9a compared directly against the SM prediction of the dimuon charge 
asymmetry. 



The first preliminary measurements of [|62|l flavour tagged B^ — t- J/^ decays at LHCb using 



36 pb ^ of data found a (statistics limited) result that agreed well with the past measurements 



of D 115911 and CDF [16011 . and was consistent with the dimuon anomaly. However a more precise 
measurement was subsequently reported at Lepton-Photon 201 1 that introduces significant tension 
into the interpretation of the pattern of data as implying new physics. This measurement reports 

Ar, ~ 0.123 ± 0.029 ± 0.008 and S^^ = 0.13 ± 0.18 ± 0.07, which gives 

' ^' ^'^ - 0.021 ± 0.024 ps-^ (70) 



1-^ 



These results are to be compared to the result derived from the updated dimuon anomaly data 
—0.20 ± 0.1 Ips^^ - demonstrating a tension at about the 2a level. 

This new LHCb result is not a direct measurement of the dimuon asymmetry, but the prospect 
for such a determination is promising. In part due to the improved mass resolution at LHCb, 
a measurement of the difference in the asymmetry in Bs mixing (a'fj and Bd mixing (aij is 
being pursued. This approach has the added benefit of detector asymmetries (but not production 
asymmetries) cancelling through taking the difference. For both a^^ and aj^, semileptonic B — )■ 
D fj,^ u decays and subsequent kaon-rich decays of D'^ — )• A'+ K^ 7r+ and D^ — )• K^ K~ 7r+ are 
used for event selection. Measuring a'^f^—a'jg also has the added benefit of producing an orthogonal 
constraint in the aj^, ai^ plane which offers strong prospects of direct confirmation or refutation 
of the dimuon anomaly in the 201 1 data set, if limitations from the production asymmetries can be 
overcome. 

At this time, available combined fits indicate a 3.3a deviation from the SM y, |21|] based on 



the older measurement and not incorporating the new LHCb and D^ measurements. We use this 
result in our analysis. In the model we have studied and attempted to link to baryogenesis, one 
expects hs = hd and a^ = a^, and the best fit values are hg = 0.255 and 2ag = 180° + 63.4° based 
on the smaller data set. 
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B. Collider prospects for 2HDMs with successful EWBG 

From a collider search perspective, the highly constrained mass spectrum and parameter space 
consistent with EWBG is a positive feature, since the predictions for collider signatures are sharp- 
ened. The typical mass scales of the new scalar states with successful baryogenesis fall into two 
categories. When we set (^ = C'b ~ in our scans we find typical values 260 < rriH ^ 360 GeV, 
350 < rriA ^ 550 GeV and 300 ^ m± < 350 GeV as is shown in fig [3l Alternatively, if we 
bias the MC search to give as large a value of C,^ = (\ as possible then we find typical masses 
275 < mj/ < 375 GeV, 475 < m^ < 575 GeV and 400 < m± < 500 GeV. In both cases, the 
simultaneous best fit value for the Higgs mass is tightly bounded around 115 ^ uih ^ 130 GeV. 
Further, in both cases A5 ~ for successful EWGB, so the new scalars, when produced, do not 
dominantly decay to pairs of Higgs particles and the mass and VEV eigenstates of the new scalars 
approximately coincide. This makes discovery of these states extremely challenging when the 
scalar masses are above the 1 1 mass threshold as we will show. SM Higgs collider phenomenology 
is not significantly affected in this scenario, and we focus our discussion on the possible collider 
signatures ofmn, ttia, m±. The challenges of discovering charged scalars at LHC are well known 
and are not alleviated in this model, so we focus our discussion on the neutral scalars. 

1. Tevatron Physics 

The lower end of the above mass intervals are well within the kinematic reach of the Tevatron 
and the simplest topology of interest at colliders is a resonant s-channel tree-level exchange of the 
new scalar states. The partonic tree-level differential cross section for mH-, mA is given by 

da _ I 1 (g-m^-rra^)(g-m|-m^) ,^c//g|2|^c//g|2 n^^ 

\Jij \ \Jkl \ ■ V'^'' 



di 192 7rs 1^2 (s_m|/^)2 + rc„,e.m 



Sr/Si ''^H/A 



We denote partonic variables with hat superscripts and the masses are labeled with 1, 2 for initial 
states and 3,4 for final states. Here the three-momenta of the initial SM states in the center of 
mass (CM) frame are p. The general coupling of the scalar fields to initial state flavours i,j and 
final state flavours A;, / is given by 

where y^ = y/2m^ /v. The expression for fki is identical. Considering parton distribution 
function (PDF) and CKM suppression, the largest effect of these exchanges at the Tevatron is 



fU/D ^ U/D 
J ij tii 



(72) 
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expected to be on the production of 1 1 final states. Specializing to this case, neglecting all initial 
state masses in the kinematic expressions, and taking yt ~ 1, we find the partonic cross section for 
i, j initial states that must be convoluted with PDF's is given by 

- _ _J}^ (-^ ~ ^"^t) I f |2 /'a0n2 /70N 

"- - 192 vr (. - </J2 + r|^/,^ m^/^ "^-" ^^* ^ " ^''^ 

It is interesting that there is an excess in the lowest mass bin between 350 and 400 GeV in the top 
quark invariant mass distribution. This range approximately coincides with the masses preferred 
for EWBG in our model when (\ ^ 0. However, we find that the maximum enhancement of 
the top quark invariant mass distribution in the relevant bin is too small by a factor of ~ 100 for 
couplings consistent with particle physics constraints and successful EWBG. The effects of tree- 
level exchanges of neutral scalars consistent with EWBG in other observables at the Tevatron are 
similarly too small to be observed with the current data set. 

2. LHC Physics 

The production of the neutral scalars through the process gg — )■ sr/j at the Tevatron and LHC is 
a promising channel for this scenario. When we bias for large values of ^^° the process pp —> sjibb 



at LHC is also promising ll25ll . Conversely, for the parameter sets consistent with EWBG we 
find that the production cross sections for pp — )■ W^ — )■ s"^ sr/i at the Tevatron and pp ^ 
W^ — )■ s^ sr/i at LHC are ^ 10"'^ fb. In Table |V] we display these cross sections for several 
2HDMs consistent with EWGB: five sets in which (\ ^ 0, and five sets in which (\ , using 
the stronger particle physics constraints that arise from assuming that / = 1 in the relation (l24l) 
between (\ and C,^. Table IVTl gives the corresponding results for models that can satisfy both 
the requirements of baryogenesis and the Tevatron dimuon anomaly, using the relaxed constraints 
from taking / = 0.5. 

For the parameter sets consistent with EWBG, we find that the production through gluon fusion 
indicates that LHC's current run (with > lfb~^ on tape) should have produced hundreds of events. 



The branching ratios of the neutral scalars are determined using the results summarized in 116311 . 
Recall our convention that the field Sdoesnot get a VEV, because of this and the fact that IA5I ^ 1, 
the mass eigenstate fields only get small VEVs of the order 

, , , v^X? , ,, v^Xi 

{sr} = vesj, = —2 -2, {sj) = ves, - 



2 2 ' \ 1 / ^^01 2 2 ■ 
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0.93 
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4.2 



TABLE V: Sample of production cross sections for models with successful EWGB, but insufficient dimuon 
excess. We use MSTW2008 PDF's in determining the cross sections and use the renormalization scale 
H = mfj/j^/2 for the PDF's. The cross sections are evaluated for LHC at an operating energy of 7TeV at 



leading order. We use the results of |[63 , 



6411 for the LO cross sections. Here we have assumed the generic 



relation (|24l) with / = 1 between ('^ and (^. 
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TABLE VI: Same as table |Vl but now with sufficient dimuon excess, allowed by mild relaxation of the 
relationship (l24l ) between C,^ and C^ from / = lto/ = 0.5. 



When EWBG occurs, these vacuum expectation values are (s'^) ^ {s'j) ~ 10 GeV. This sig- 
nificantly suppresses decays to WW, ZZ making these scalars difficult to discover at LHC. The 
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largest decay widths are typically to fermions with 

r{sR/j^ff) 



GpNc 2|^0|2 on 



(74) 



8^2 vr 

and /3 = (1 — 4mj/m^/^)^/^. Here n = 3 for m// and n = 1 for m^i. The neutral scalars decay 
almost exclusively into 1 1 final states when they are above the 1 1 threshold. Other phenomenolog- 
ically interesting final states due to their clean collider signatures are 77, ZZ, WW and r+ r~. 
We show in Fig. (fT2l) the dominant branching ratios as a function of the scalar masses for typi- 
cal parameter values. We also show the branching ratios to EW final states that are potentially 
promising for discovery. 
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FIG. 12: Branching ratios for the heavy neutral scalars. The branching is almost exclusively to VF VF and 
Z Z below the top quark threshold and to 1 1 above the top threshold. Here we have neglected flavour- 
changing decays and decays to hh, the later being highly suppressed compared to the leading decay to 
W W and Z Z or tt and challenging to discover. Each curve is labeled by the final state particle X for the 
decay s^// — )• XX. We have taken (^ = 0.5, mt = 173.1 and (^ = 2m'^/v'^. Flavour changing decays to 
t c or tu require an insertion of 2ml /v'^ ^^'^ further CKM suppression. They are given essentially by CKM 
suppressed rescalings of the shown Br (6 6). 

Below the t i threshold, strong discovery prospects exist through decays to VF ly and Z Z 
using standard searches. Typically in these scenarios at least one heavy neutral scalar is below 
the t i threshold. However above the t i threshold, discovery prospects are much dimmer. When 
TT^H/A > 2 mj large numbers of resonant t i pairs are produced. Unfortunately the t i production 
from QCD is far larger because the signal process is loop-suppressed. For example, for sj —^ ti 
with a mass scale of ttia ~ 500 GeV, the approximate cross section for cr{g g — )• s/) is 0.2 fb 



48 

while the QCD 1 1 background with Mti ^ 500 GeV is given by ~ 10 pb at leading order. The first 



measurements of the 1 1 invariant mass spectrum at LHC have recently appeared 116511 . However it 
is unlikely that the resonant production of 1 1 due to these heavy scalars can be uncovered in the 
tt invariant mass spectrum for quite some time. Precision studies would have to be performed, 
possibly using the angular distributions of the reconstructed 1 1 pair since this signal is through a 
scalar resonance. After event selection cuts, the reconstructed partonic t i spectra has only 5 events 
in the Mti ^ 500 GeV mass region in this study. Thus total event selection cut suppression of the 
background and signal is ~ 10^^. 

Due to the difficulty of overcoming the 1 1 QCD background, and the lack of efficient light 
flavour tagging at LHC to uncover, for example, CKM suppressed t c decays, utilizing a highly 
suppressed branching ratio such as to r+ r^ might be considered more hopeful. The dominant 
background'"* in this case is the inclusive process pp — )■ Z X — )• t'^t^ X, whose tree-level 
electroweak contribution has an antiquark initial-state PDF suppression at LHC. The partonic cross 
section for this background is given by 



J\2 , (J\'2 



^^ 48 TT sin^ ^VF cos^ ^VF {s - mlf + mlVl 



(75) 



Here gv,gA are the usual electroweak couplings that depend on the isospin of the initial state 
fermions /. We use Tz = 2.5 GeV and m^ = 9L2GeV. We use hat superscripts to denote 
partonic quantities. The expression we use to implement cuts on the 2 — ?> 2 processes of interest, 
in fb units, is given by 

da 0.3894 X 10^^ M., fy^ p 

■77^ = / dye dz^ff{Rey^)fj{Re y^)af (76) 

UMt-t Shad J-yO^ J-zo f 

where R = ^J M'^^jshad, fj, fj are the 2008 MSTW parton densities of the initial state quarks. 
Shad is the hadronic center of mass energy, and Yb is the boost rapidity of the subprocess frame. 
We treat the highly boosted r's as r jets required to have rapidity |yi_2 1 < Y^ax with cuts on px of 
p??^" so that 

1/°, = Min[i;„,,-Log[M2^/s]/2], 

z° = Min[y/l - 4 {pf^Y/M^^^, tanh(i;„,, - ly^l)]. (77) 



'^ This is most difficult to significantly reduce with selection cuts. A background of tt also exists with r pairs 
produced in the decay of the tops. The latter should be more reducible due to its associated missing energy and 
multi-jet topology. 
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This result for hadronic cuts on 2 — )> 2 processes is consistent with [|66n . We find that the d a/dM^-r 
invariant mass spectrum has a peak of resonantly produced r+ r^ final states typically 10~^ the 
background spectrum when apr cut of 10 GeV is applied to the leptons and the r jets are restricted 
to have a rapidity of < 3. Thus discovery in this channel is unlikely. Above the ti threshold, 
relying on decays to VT TV and Z Z offers some prospects for discovery. 

VIII. DISCUSSION AND CONCLUSIONS 

In this paper we have given an exhaustive account of EWBG in two-Higgs doublet models. 
Our study, while focusing on MFV Yukawa couplings and a simultaneous explanation of the D^ 
dimuon excess, is general and encompasses earlier popular choices of Yukawa couplings, such 
as type I and II models. The MFV two Higgs doublet model also offers a controlled theoretical 
framework where new CP violating phases can appear at the EW scale without introducing in- 
consistencies with a multitude of other flavour measurements. Unlike previous works, we have 
considered the most general set of couplings in the Higgs potential and have made no simplifying 
assumptions to reduce the parameter space. Instead we used a Markov Chain Monte Carlo to ex- 
plore the full space of models, aiming toward those with a strong first order phase transition and 
large CP violation in the bubble walls. 

Despite this, we find it more difficult to obtain a large enough baryon asymmetry than previous 
studies indicated. This is partly because we have imposed more phenomenological constraints; 
in addition to accelerator mass bounds and the neutron EDM that were considered in previous 
studies, we included EWPD, Rf„ b ^ s'j and perturbativity of couplings up to a minimum Landau 
pole scale of 1 — 2 TeV In addition, we have avoided simplifying assumptions about the shape 
of the bubble wall profile that artificially increased the baryon asymmetry in previous work. For 
every model in the chain from the Monte Carlo, we numerically determined the actual solutions for 
the bubble wall profiles and the CP-violating phase 9{z) that sources baryogenesis. In the process 
we elucidated a dynamical mechanism by which 9(z) tends to be suppressed: the existence of 
nonzero 9{z) depends upon exciting the field S whose VEV is zero at T = 0. Thus there is 
finite-temperature suppression factor inherent in 9{z) which would be overlooked in a seemingly 
reasonable guess for its profile. We numerically solved the Boltzmann transport equations for each 
profile to determine the baryon asymmetry. 

The effective potential approach we have employed has theoretical uncertainties that are diffi- 
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cult to quantify without going through the arduous step of adding two-loop corrections. Ultimately 
a lattice study would be welcome for verifying the details of the phase transition properties, but 
this is an even more ambitious undertaking.'^ A more tractable improvement on our treatment 
would be to rederive the transport equations and related coefficients, which we have taken from 
ref. i2m without modification. In particular, it is possible that including extra species in the net- 
work of equations could enhance the baryon asymmetry. Tau leptons, although weakly coupled, 
diffuse farther than other species, and could thus carry the CP asymmetry deeper into the sym- 
metric phase and boost the integral (l66l) determining n^.'^ A crucial issue to resolve concerns the 



correct form of the source term in the transport equations. Ref. yJO has recently claimed that 
2HDMs can more readily produce a large enough baryon asymmetry (as well as address B meson 
anomalies including the one we have considered) than we find. This is in part due to the use of 
a source term that is first order in derivatives, rather than the semiclassical force term with two 
derivatives that we have adopted. 

Given these theoretical uncertainties in computing the baryon asymmetry, it seems reasonable 
to believe that the true prediction for a given model could differ by a factor of 2 from our de- 
terminations, despite the pains we have taken to get accurate results. At the same time, there 
is uncertainty in the value of the new Yukawa coupling parameter 77^ needed to match the D^ 
anomaly; from table 17 of [4], we can estimate that the 1-cr error band in hg of eq. (fT6l) is 



hq = 0.258 ^QgQ using the method advocated in [|67ll . Since hq ~ ('7d)^, this translates to the 
relative error 5?7^/|77^| = ^q |^54- From these considerations we can give a better picture of the 
tension between explaining the dimuon excess and baryogenesis as shown in fig. \T3\ It will be 
important to reduce the systematic errors in the prediction of the baryon asymmetry to know how 
serious the apparent tension really is. 

One advantage of having a small allowed region of parameter space is that we can make some- 
what distinctive predictions for collider searches. Typically the discovered mass spectrum is a 
light Higgs, and one or more heavy scalars below the t i mass threshold with branching ratios as 
shown in fig. [I2l LHC can discover such a spectrum in the near future. This scenario for EWBG 
is falsifiable for most of the remaining allowed parameter space. 



'^ Ref. 1681] has also promoted a method to extract Vc/Tc from the effective potential in a gauge-invariant manner. 
'^ We thank S. Tulin for pointing out this possibility 
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FIG. 13: Scatter plot of \rj'^\ normalized to the observed value (as needed to explain the Y)(j) dimuon excess) 
versus the baryon asymmetry relative to its observed value. The 1, 2 and 3-o- eiTor ellipses as described in 
the text are shown. Black points are for / = 1, magenta for / = 0.5, green for / = 0.2, referring to / 
defined in eq. (l24l) . 
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Appendix A: Counterterms 

Denoting derivatives of the Coleman- Weinberg potential Vcw^ evaluated at the tree-level VEVs, 
by Vi^^ and Vi^^y, 

(Al) 

(A2) 
(A3) 
(A4) 
(A5) 
(A6) 
(A7) 
(A8) 
(A9) 

As described in the text, the Goldstone boson contributions to Vcw must be IR-regulated in order 
to give meaningful results. 

Appendix B: Field-dependent masses 

The neutral components of the complex scalar doublets can be written as 

H = -={h + zhi), S = ^{sR + isi) (Bl) 

In our evaluation of the effective potential we are free to make a global SU(2)i transformation 
that sets hj = so that V^cfj depends only upon h, sji,sj. The resulting 4x4 neutral mass matrix, 



5Xv^ = 


U-Via-V^m) 


5^^ = 


U-Vi^-v^m) 


6Xfv^ = 


2i~^l,SRSR + yi,SlSj) 


6Xiv^ = 


^l,SBSi 


5ml = 


-W^i'^R^R + ^1,^/^/) 


5m\p, = 


2i^l,hsR - -Vi^sr) 


6mlj = 


U-VlMr + !^mJ 


6X^v^ = 


-y^hsR + ^v^i,sfl 


5XW = 


-Vi^hsj + lVl,sii 
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including finite-T Debye mass corrections, has components given by 

VH,h = \ {Sh' - v') + K^i + ^2)(4 + sj) + A3 (4 - sj) + 3h{Xfsn + Xfsj) 

+ (3 A + 4Ai + 2A2 + 6y^ + ^g' + y") g (B2) 

Vh,,h, = ^{h'- v') + i(Ai + X,){s% + sj) - A3 (4 - ^?) + MAf^iJ + Af .,) 

+ (3A + 4Ai + 2 A2 + 6y2 + 1^2 ^ 1^/2) g (33) 

Vsu,sn = m^ + 1 (Ai + A2 + 2A3) h^ + 3Af hsR - Xi hsj + Aq (34 + s]) 

+ (4Ai + 2A2 + 12A6 + 6|r/i|2 + |(72 + |^'2)|f ^34) 

F,,,,, = m^ + I (Ai + A2 - 2A3) h^ - 3X{ hsi + Af /iSR + As (3s? + 4) 

+ (4Ai + 2A2 + 12A6 + 6|r/i|2 + ^g^ + y^)^ (B5) 

V^,h, = X^iSRSi + h{X^si-XisR) (B6) 

V?^,.« = m^^ + (Ai + A2 + 2A3) hsR + |Af (3s^ + s^) - X{sr si + l X^h^ 

+ 6{Xf + Xf + yt\vt\ cos ipr,)^ (B7) 

Vh,sj = -mil + (Ai + A2 - 2X3)hsi - ^Xi (4 + 3s?) + XfsR s/ + f A^ h^ 

+ Q{Xi-Xi- yt\vt\ sin ^r,)^ (B8) 

Vhi,s^ = mlj + 2 A3 /is, + ^Xi (34 + s?) + XfsR sj-^Xi K" (B9) 

K,,s, = m^^ + 2 A3 /isfi + |Af (4 + 3s?) + A^ Sij s/ + I Af /i" 

+ 6(Af + Af + 2/t|r/i|cos(^^)g (BIO) 

Kh,s, = KX'^si -XisR) + 2XQSiSR (BU) 

Here rjt = \'r]t\e^'^'^ and we have denoted the real and imaginary parts of complex masses and 
couplings by ml = mlj^ + imlj, etc. The mass matrix must be numerically diagonalized to find 
the mass squared eigenvalues. 
The charged scalar mass matrix elements are given by 

Vh+,h- = ^{h'-v') + lX^{sl + sj) + h{XfsR + Xisj) 

+ (3A + 4Ai + 2A2 + 6y^ + |/ + Ig")^ (B12) 

Vh+,s- = (ml)* + |A2/i(sR + isi) + XshisR - tsj) + ^XHsl + s?) + ^X^h^ 

+ 6(A5 + A: + 2/ir/;)g, (B13) 
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Vs+,h- = K+,s- (B14) 

K+,s- = ml + iAi/i^ + /i(Af sk - \{si) + \{s\ - s]) 

+ (4Ai + 2A2 + 12A6 + 6|r/t|2 + |/ + f^'^^g ^ (515) 



and the charged scalar mass eigenvalues are 

1 



^l = i^{yh+,h- + K+,s-) ± l^^J{Vh+,h~ - Vs+,s-Y + 4|V;+,,-|2 . (B16) 



The gauge boson masses in the basis of W^, B"' are given by 



ml,„,..,, = {lg^(f)' + 2g^T^6^,\i)6,j6^0 (B17) 



W°'W^ 



m 



2 

W°'l 



= {l99'<P')Sap (B18) 

rnlo^B^ = (i^V + 2^"T2(5,,|,)5„^ (B19) 

Here 0^ = /i^ + s|. + sj denotes the total electroweak-breaking VEV and 5a,\\ indicates that only 
the longitudinal polarizations get a Debye correction at this order. 

Finally, the top quark mass is given by 

^t = I {{yth + \r]t\{s R cos Lfn - SjsmLpn)f + \r]t\'^ (sr sin Lpn + Sj cos LPnY) . (B20) 



m. 



Appendix C: Beta functions 



The scaled beta functions, generalizing the results of ref. [|32D, are as follows. We have defined 
A^'^ to be the real or imaginary part, respectively, of A„. By "scaled," we mean that (3i = /3j/167r^ 
in terms of the conventionally normalized function, (3x = dX/d In /i, where jj, is the renormalization 
scale. 

4 = 6 A^ + 8A2 + 8A1 A2 + AXl + 16A2 + 48|A5p + ^ (3/ + g'^ + 2g^ g'^) 

-3A(3/ + (7'2-%2)_l2|/,^ 
/3ai = (A + 4 As) (3Ai + A2) + 4A? + 2A^ + SA^ + 4IA5I' + 16 (Af Af + A^ A^) + 4|A4|', 
+\ i^g' + g" - 2g'g'') - 3Ai (3^^ + g'2 _ 3^2 _ 2|^^|2 _ 2|^^|2) _ 12^2 j^^p^ 

/3a2 = (A + 4 Ae) A2 + 8 Ai A2 + iXJ + 32 A^ + IOIA5I' + 4 (Af Af + A^ A^) + IOIA4I', 
+3^V' - 3A2 {3g^ + g'^ - 2yl - 2W\^ - 2\i]^f) - \2yl\^t? , 
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^3 = (A + 4A6)A3 + 8AiA3 + 12A2A3 + 5|A5p + 2(AfAf + A^A3+5|A4p, 

-3A3 {3g^ + 9" - 2yl - 2|r/,|2 - 2\t^,\^) - Qy^M\\ 
/3j = 24A6 Af + 6Ai(Af + Af ) + 4A2 Af + 8A2 Af + 4A3 Af + 2OA3 Af 

-3 Af (3/ + g'' - vl - ZW\^ - ZW\^) - 12 {(r^ff + {r,\f (r^f )) y„ 
/3{^ = 24A6 \i + 6Ai(-A^ + A^) - 4A2 A^ + 8A2 A^ - 4A3 A^ - 2OA3 A^ 

-3 \{ (3/ + g'^ - y\ - 3|ry,p - ?.W\^) - 12 ^f + (ryf )^ (ryf )) y„ 
/3j = 6A Af + 6Ai(Af + Af ) + 8A2 Af + 4A2 Af + 2OA3 Af + 4A3 Af 

-3 Af (3/ + (?'2 - 3y2 - |ry,|2 _ |^,|2) _ llv^fyl 
hi, = 6XXI + 6Ai(A^ - A^) + 8A2 A^ - 4A2 A^ - 2OA3 A^ + 4A3 A^ 

-3 A^ (3/ + g'^ - 3y^ - \7^t\' - \Vb\') - I'^vlvl 
/3a6 = 24A2 + 2\\ + 2Ai A2 + A^ + AXl + 12A2 + ^ {?,g^ + g'^ + 2g^g'^) 

o 

-3A6 (3/ + g'^ - 4|r/i|2 - A\r],\'') - 6|r/i|^ - Q\r],\\ (CI) 
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